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Abstract 

Phylogenetic relationships of snailftshes of the family Liparidae were analyzed on the basis of two sets of molecular 
sequence data: one from the mitochondrial DNA cytochrome c oxidase subunit one gene (COI) and another from 
restriction-site associated genome-wide sequences (RADseq). The analysis of COI sequence data from at least 122 species 
of 18 genera from the Pacific, Atlantic, and Southern oceans resulted in a moderately well-resolved phylogeny among 
the major clades, albeit with significant polytomy among central clades. Nectoliparis was the sister of all other members 
of the family, followed by Liparis. Liparis, Careproctus, and Pamliparis were paraphyletic. Liparis was recovered in 
two closely related clades, with L. fucensis sister of all other liparids except Nectoliparis, and both Careproctus and 
Pamliparis were each recovered among at least three widely separated clades. The RADseq analysis of 26 species of 11 
genera from the eastern North Pacific strongly confirmed the overall results of the COI analysis, with the exception of the 
paraphyly of Liparis due to the absence of L. fucensis. Our results show that the pelvic disc has been independently lost 
multiple times and the pectoral-fin girdle has been independently reduced in multiple lineages. 

Key words: Teleostei, molecular, morphology, Liparis, Careproctus, Pamliparis, phylogeny 


Introduction 

Snailftshes of the family Liparidae are a well-defined group of diverse cottoid fishes characterized most promi¬ 
nently by their scale-less, typically naked, elongate bodies, and by the development of a sucking disk fonned by 
highly modified pelvic fins. The family is bipolar in distribution and is one of the four most-rapidly speciating 
clades of marine fishes, all found in high-latitudes (Rabosky et al. 2018). The Liparidae comprises over 430 species 
(Eschmeyer & Fong 2017) allocated to about 32 genera (Chernova et al. 2004; Orr 2004; Balushkin 2012; Stein 
2015; Eschmeyer et al. 2017). Two of its most diverse genera, Careproctus and Pamliparis, are recognized in both 
the Northern Hemisphere (northern species) and the Southern Hemisphere (southern species), while a third, Liparis, 
is restricted to the north. Presently, 214 species are known from the north, primarily the North Pacific, whereas 262 
species are recorded from the south, primarily the Southern Ocean, and many new species continue to be described 
both north and south (e.g., Chernova 2005ab, 2014; Stein 2006, 2012, 2016; Orr 2102, 2016; Eschmeyer & Fong, 
2017). Liparids are particularly diverse in North Pacific waters, especially off Alaska where over 85 species have 
either been described or are known but as yet undescribed (Mecklenburg et al. 2002; Orr 2004; Orr & Busby 2006; 
Orr & Maslenikov 2007; Orr 2012, 2016; Orr et al. 2015; Orr, unpublished). They are abundant as well as diverse, 
as documented in trawl surveys conducted throughout a wide range of depths, especially in surveys of the Aleutian 
Islands (Raring et al. 2016) and upper continental slope of the eastern Bering Sea (Hoff 2016). 

While bipolar in distribution, and essentially limited to cool-temperate and cold waters, liparids are widely 
distributed from the shallowest waters of the intertidal region to near the deepest parts of the ocean. The recently 
described Mariana Snailfish, Pseudoliparis swirei, has been documented at depths to 8178 m and is now recognized 
as the deepest dwelling vertebrate thus far confirmed (Linley et al. 2016; Gerringer et al. 2017a; JAMSTEC 2017). 
In addition, although not yet captured and thus not formally described, the Ethereal Snailfish was filmed in the 
Mariana Trench at a depth of 8,145 m, only 33 m shallower (Linley et al. 2016; Gerringer et al. 2017a). Most spe¬ 
cies are demersal or benthic, but many species are found also in the midwater. Their precise contributions to marine 
ecosystems are unclear, but liparids are important as prey and as predators in many ecosystems (Bansode et al. 2014; 
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Gerringer et al. 2017b; Whitehouse et al. 2017). Among their unusual behavioral characteristics, species of Care- 
proctus are known to lay their eggs within the carapace of lithodid king crabs (Stein 1980; Gardner et al. 2016). 

Phylogenetic status. Historically placed within the Scorpaeniformes among the sculpins (Nelson et al. 2016), 
the Liparidae was recently reallocated to the perciform suborder Cottoidei by Smith & Busby (2014; = infraorder 
Cottales of Betancur-R et al. 2017). Its close relationship with the Cyclopteridae, as originally proposed by Peter 
Artedi (1738; see also Walbaum, 1792), followed later by Cuvier (1817) and Garman (1892), has been unchal¬ 
lenged, and both families are unquestionably monophyletic (Near et al. 2012; Smith & Busby 2014). While sys¬ 
tematic progress has been made with the Cyclopteridae (Kai et al. 2015; Oku et al. 2017; Stevenson et al. 2017), 
the phylogenetic relationships of the far more diverse Liparidae are poorly understood (Smith & Busby 2014) and 
the monophyly of major genera within the family has been questioned (Andriashev 1986; Rnudsen & Moller 2007; 
Balushkin 2012). Previous studies have been limited primarily by inadequate taxon sampling of this morphologi¬ 
cally diverse and species-rich group, as well as by the challenges of working with the morphology of these poorly 
ossified, reductive, and easily damaged fishes. 

Early works focused on the few species known at the time and gave only general ideas of their evolution and 
relationships: Garman (1892) provided the first detailed anatomical descriptions and reviewed all 18 known spe¬ 
cies in the three genera— Liparis, Careproctus, and Paraliparis —recognized at the time; Goode & Bean (1896) 
reviewed 10 deep-sea species in five genera; and Jordan & Evermann (1898) included all 39 species known among 
10 genera in North American waters. Burke’s (1930) monumental revisionary work on the family summarized all 
that was known at the time about the family, its 114 species in 13 genera, and its evolution. As in the works of pre¬ 
vious authors, Burke described the family as being closely related to the Cyclopteridae and more distantly to the 
broadly defined Cottidae, supporting the early consensus about evolution within the group: that is, their evolution 
from shallow-water, large-disc species of Liparis, through disc reduction among deeper dwelling Careproctus, to 
the complete loss of the pelvic disc among deeper dwelling and semi-pelagic species of Paraliparis and related 
genera. Thus, genera were arranged in phylogenetic order from least to most derived, with those more basal bearing 
a disc (i.e., Liparis and Polypera, Careproctus and Temnocora, Ciystallichthvs and Crystallias, and Gyrinichthys ) 
and those more derived lacking or having a highly reduced disc (i.e., Elassodiscus, Paraliparis, Rhinoliparis, Acan- 
tholiparis, Nectoliparis, and Lipariscus ). 

Overall, later analyses supported these general patterns, but some notable exceptions were demonstrated. Es¬ 
pecially as the full diversity of the Liparidae has become more clearly understood, wide gaps in our understanding 
of the phylogenetics have been revealed. Several modern phylogenetic analyses have been published based on 
morphological and molecular data (Figs. 1—7). As in other complex groups, analyses have moved from those using 
strictly morphological data, to combined analyses including both molecular and morphological data, and most re¬ 
cently to strictly molecular data. Taxon and individual sampling have generally increased among species and genera 
with each analysis. 

The earliest analyses used morphological characters exclusively. In the first modern phylogenetic analysis, 
Kido (1988; Fig. 1) examined 60 species in about nine genera using osteology and external morphology, focusing 
on species of Japanese waters augmented by selected species from the North Atlantic. In challenging aspects of 
Kido’s (1988) results, particularly the position of Nectoliparis as sister of all other liparids, Balushkin (1996; Fig. 2) 
attempted to describe the relationship of his new monotypic genus, Palmoliparis, among the other 26 known liparid 
genera. While preliminary and supported by only seven osteological and external characters, his analysis was the 
first attempt to place all known genera within a cladistic scheme based on well-recognized character transforma¬ 
tions and current interpretations of liparid evolution (Andriashev 1998; Andriashev & Stein 1998). In both of these 
morphological analyses, Liparis was monophyletic and Careproctus was paraphyletic (in the reanalysis of Kido’s 
1988 data by Knudsen et al. 2007). Paraliparis was also paraphyletic in Kido’s (1988) analysis; its status was not 
examined by Balushkin (1996). 

The first phylogenetic analysis to use molecular data was that of Knudsen et al. (2007; Fig. 3), who conducted a 
combined analysis of two mitochondrial DNA (mtDNA) datasets of over 1000 bp with morphological data, based on 
24 species in nine North Atlantic and North Pacific genera. With this limited taxon sampling, they found all genera 
except Paraliparis to be monophyletic. 

More recent species-level analyses have used only mitochondrial sequence data, primarily the cytochrome c 
oxidase subunit I gene (COI), and have focused on the identification of species. As collections of tissues have in¬ 
creased, numbers of species, genera, and individuals have grown with each analysis. 

Steinke et al. (2009; Fig. 4) analyzed sequence divergences in a 650 bp alignment of COI of 78 tissues of 19 
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species of eight Canadian Pacific genera. Their goal was to evaluate the utility of COI for distinguishing species, 
although misidentifications hampered the accuracy of their species-level conclusions. In their neighbor-joining 
tree, Liparis was monophyletic and both Paraliparis and Careproctus were paraphyletic. Similarly, Duhamel et al. 
(2010; Fig. 5) focused on Antarctic species and significantly increased the overall number of species and geographic 
scope sampled by analyzing 157 tissues of 46 species in nine genera. With much new material of Antarctic species 
and 103 COI sequences 668 bp in length from other geographic areas taken from online databases, their Bayesian 
analysis found not only Paraliparis and Careproctus to be paraphyletic but, for the first time, Liparis as well. 

- Paraliparis pectoralis 

- Paraliparis grandis 

- Paraliparis dactylosus 

- Paraliparis ulochir 

- Rhinoliparis barbulifer 

- Rhinoliparis attenuatus 

- Paraliparis opercularis 

_ _ Paraliparis rosaceus 

- Rhodichthys regina 

- - - Paraliparis garmani 

- Paraliparis bathybius 

- Lipariscus nanus 

- Elassodiscus tremebundus 

- Careproctus group 

- Careproctus mederi group 

- Careproctus reinhardti 

- Careproctus dentatus 

- Osteodiscus 

- Crystallichthys cyclospilus 

- Crystallichthys mirabilis 

- Crystallichthys matsushimae 

- Liparis montagui group 

- Liparis pulchellus 

- Liparis group 

- Liparis fabricii 

- Liparis greeni 

- Nectoliparis pelagicus 

FIGURE 1 . Phylogenetic hypothesis of Kido (1988, after fig. 20), derived from a maximum parsimony analysis of morphologi¬ 
cal data, using 34 osteological and external characters, for 60 liparid species. 
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- Nectoliparis 

- Lipariscus 

- Psednos 

- - Edentoliparis 

- Odontoliparis 

- Genioliparis 

- Acantholiparis 

- Rhinoliparis 

- Rhodichthys 

- Paraliparis 

- Elassodiscus 

- Careproctus 

- Osteodiscus 

- Notoliparis 

- Pseudoliparis 

- Pseudonotoliparis 

- Squaloliparis 

- Eknomoliparis 

- Allocareproctus 

- Temnocora 

- Gyrinichthys 

- Careproctus 

- Careproctus 

- Palmoliparis 

- Crystallias 

- Crystallichthys 

- Polypera 

- Liparis 

FIGURE 2. Phylogenetic hypothesis of Balushkin (1996, after fig. 4), derived from a manual cladistic analysis of morphologi¬ 
cal data, including seven osteological and external characters, for 26 liparid genera. 
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Psednos micruroides 


1.00 


1.00 

1.00 

1.00 






Psednos melanocephalus 
Psednos groenlandicus 
Rhodichthys regina 


1.00 


0.97 


0.87 


0.98 


1.00 


1.00 


0.74 


1.00 


1.00 


1.00 


Paraliparis garmani 
Paraliparis sp. 1 

Paraliparis bathybius 
Paraliparis copei 

Rhinoliparis barbulifer 
Elassodiscus caudatus 

Elassodiscus tremebundus 
Careproctus cypselurus 

Careproctus melanurus 
Careproctus kidoi 

Careproctus rastrinus 
Careproctus reinhardti 

Crystallichthys cyclospilus 
Liparis montagui 

Liparis liparis 
Liparis inquilinus 

Liparis fabricii 
Liparis gibbus 
Liparis pulchellus 
Nectoliparis pelagicus 


0.98 


0.92 


Eumicrotremus derjugini 
Eumicrotremus spinosus 

Cyclopterus lumpus 
Aptocyclus ventricosus 


FIGURE 3. Majority-rule (50%) consensus tree of Knudsen et al. (2007, after fig. 3), derived from a Bayesian analysis of three 
combined datasets composed of mitochondrial DNA (16S and cytochrome b ) and morphological data for 24 liparid species. Tree 
is rooted with species of the Cyclopteridae. Posterior probabilities are above branches. 
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Paraliparis sp. (P. sp. cf. rosaceus) 


83 


89 


98 


99 


Paraliparis paucidens 


100 


- Lipariscus nanus 

Paraliparis melanobranchius (P. megalopus) 

Paraliparis dactylosus (P. ulochir) 

- Paraliparis rosaceus 


Rhinoliparis attenuatus 

Paraliparis pectoralis I 

100 


— Paraliparis pectoralis 2 
Careproctus georgianus (C. ovigerus) 
- Acantholiparis opercularis 


Careproctus attenuatus (C. longifilis) 


100 


100 


— Careproctus cypselurus 1 (Osteodiscus cascadiae ) 

— Careproctus cypselurus 2 (Osteodiscus cascadiae) 
Elassodiscus caudatus (Paraliparis dactylosus) 

- Careproctus canus (Careproctus n. sp. 

- Careproctus melanurus 


97 


- Careproctus furcellus (C. cypselurus) 
Liparis dennyi 

- Liparis pulchellus 

- Nectoliparis pelagicus 


0.02 


FIGURE 4. Unrooted neighbor-joining tree of Steinke et al. (2009, after fig. 4), derived from cluster analysis of a 650 bp align¬ 
ment of cytochrome c oxidase subunit 1 gene (COI) sequences for 78 samples of 19 liparid species. Bootstrap values >80 are 
above branches leading to multiple species. Corrected identifications based on our study are in parentheses. 
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Careproctus georgianus (C. ovigerus) 

— Careproctus crozetensis 

Elassodiscus caudatus (Paraliparis dactylosus) 

- Careproctus canus (Careproctus n. sp.J 

_^|— Careproctus georgianus 


0. 601 - 

099 I— 


Paraliparis sp. 

- Careproctus melanurus 

Careproctus furcellus (C. cypselurus) 


Careproctus furcellus (C. cypselurus) 

— Acantholiparis opercularis 
Careproctus attenuatus (C. longifilis) 

— Careproctus cypselurus (Osteodiscus cascadiae) 

Careproctus longipectoralis 
Careproctus discoveryae 
Careproctus continentaiis 


-HI 


fC 


— Paraliparis pectoralis 
Paraliparis thaiassobathyaiis 
Careproctus georgianus 
Lipariscus nanus 


Paraliparis paucidens 
Paraliparis sp. 


Paraliparis mawsoni 


— Careproctus georgianus 

Paraliparis rosaceus 


_j_|- Pare 

L Paraliparis neelovi 


- Paraliparis bathybius 

Paraliparis dactylosus (P. ulochir) 

Paraliparis meianobranchius (P. megalopus) 
- Paraliparis wotffi 

Edentoliparis terraenovae 
— Paraliparis operculosus 
Paraliparis antarcticus 


ft 


Paraliparis sp. 


— Paraliparis leobergi 


l— Paraliparis charcoti 

- Paraliparis valentinae 

— Rhinoliparis attenuatus 


Liparis pulchellus 
- Liparis mucosus 


- Liparis rutteri 

i i— Liparis mucosus (L. dennyi) 

L Liparis dennyi 

— Liparis gibbus (L. bathyarcticus) 
Liparis fabricii 
■ Liparis florae (L. greeni) 


Liparis fucensis 


Nectoliparis pelagicus 
— Cyctopterus lumpus 


Zoarcidae 


FIGURE 5. Consensus phylogenetic tree of Duhamel etal. (2010, after fig. 3), derived from Bayesian and maximum parsimony 
analyses of a 668 bp alignment of cytochrome c oxidase subunit 1 gene (COI) sequences for 157 samples of 46 liparid species. 
Bayesian posterior probabilities are above branches that lead to multiple species. Tree is rooted with species of the Cyclopteri- 
dae and Zoarcidae. Corrected identifications based on our study are in parentheses. 
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1.0 I— 



98 
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67 
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Careproctus melanurus (C. sp. cf. melanurus) 

Careproctus cypselurus 

Careproctus furcellus 

Careproctus colletti 

Careproctus simus 

Careproctus scottae 
Careproctus phasma 
Careproctus gilberti 
Careproctus ostentum 
Careproctus comus 
Careproctus faunus 
Careproctus bowersianus 
Acantholiparis opercularis 
Paraliparis dactylosus 
Paraliparis dactylosus (P. grandis) 

Paraliparis pectoralis 
Lipariscus nanus 
Crystallichthys cyclospilus 
Rhinoliparis attenuatus 
Elassodiscus caudatus 

Elassodiscus tremebundus 
Nectoliparis pelagicus 
Liparis gibbus 


FIGURE 6. Majority rule (50%) consensus phylogenetic tree of Gardner et al. (2016, after fig. 4), derived from Bayesian infer¬ 
ence and maximum parsimony analysis of a 492 bp alignment of cytochrome c oxidase subunit 1 gene (COI) sequences of 492 
bp for 128 samples of 23 liparid species. Bootstrap values are above and Bayesian posterior probabilities are below branches that 
lead to multiple species. Tree is rooted with Liparis gibbus. Corrected identifications based on our study are in parentheses. 


Primarily concerned with identifying eggs of snailfishes found within lithodid crabs, Gardner et al. (2016; Fig. 
6) conducted a Bayesian analysis of 429 bp of COI data for 23 species and nine genera of 53 liparid adults from the 
Bering Sea, Aleutian Islands, and eastern North Pacific Ocean. In their analysis, multiple clades of Careproctus and 
Paraliparis, as well as other genera, were placed in an unresolved polytomy of their consensus tree. 

Most recently, Shen et al. (2017; Fig. 7) published the full mitochondrial genome of an undescribed deep-water 
liparid and attempted to place the species phylogenetically through the analysis of mtDNA data, including COI, 
taken from online databases. In analyses with the 16S and Cytb mtDNA data, their overall results were unsurpris¬ 
ingly similar to those of Knudsen et al. (2007), the primary source of their data. With the COI data alone for 440 bp 
of 83 species in nine genera, they produced a tree that resolved all three major genera as paraphyletic. 

Methodological approach. Despite these many previous efforts, the phylogenetic relationships of liparids 
remain poorly understood and increasing evidence reveals incongruence between liparid taxonomy and its phylo¬ 
genetic history. The explosion of readily available molecular data, especially the wide use of COI sequence data 
stored in public databases, and the development of new next-generation DNA sequencing techniques, which incor¬ 
porate more data than typically employed by Sanger sequencing techniques and encompass a greater range of sites 
throughout the genome, have provided new resources and tools for phylogenetic research (e.g., Leache et al. 2015). 
The time is thus ripe for more extensive molecular phylogenetic analyses. 

The mitochondrial COI gene is relatively easy to sample and because of widespread efforts to routinely sample 
species to capture biodiversity information, copious data are available from public databases, such as the Barcode of 
Life Database (BOLD; Ratnasingham & Hebert 2007) and GenBank (Benson et al. 2005). However, the gene may 
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have insufficient phylogenetic signal to resolve deeper evolutionary divergences in many groups (Hajibabaei et al. 
2006, 2007). Coupled with limited taxon sampling, it may also be positively misleading because of homoplasy and 
long branch lengths among diverse taxa (Felsenstein 1978; Anderson & Swofford 2004). 

To overcome some of the limitations of using COI alone, we complemented our analysis with data derived 
from restriction-site associated DNA sequencing, or RADseq (Baird et al. 2008), which describes a range of related 
techniques that use restriction enzymes to target for sequencing a subset of genomic regions flanking the cut sites. 
The target regions are barcode-tagged, isolated, pooled, and amplified, and next-generation sequencing techniques 
are used to generate DNA sequence data from flanking regions adjacent to restriction enzyme cut sites. Targeted 
regions include coding and non-coding loci from throughout the genome and, because restriction cut-sites are gen¬ 
erally conserved, individuals of the same or closely related species typically share most loci (Andrews et al. 2016). 
Although taxon sampling may be limited because more intensive effort and higher quality DNA is required than in 
a single-gene approach, RADseq utilizes hundreds more loci than could previously be surveyed in single- or multi¬ 
gene approaches and generates a wealth of useful data for phylogenomics. It therefore has the potential to create 
more fully resolved and better supported phylogenetic hypotheses, as well as better estimates of historical diver¬ 
gence times, than are possible with traditional studies (Eaton & Ree 2013; Cruaud et al. 2014; Longo & Bernardi 
2015; Tariel et al. 2016). 

RADseq has been successfully applied to a wide variety of species and for a variety of uses. A single RADseq 
library can be used for a wide range of applications by filtering the data at appropriate levels. Like whole-genome 
sequencing, its uses range from genome mapping (Baxter et al. 2011) to population identification (Davey & Blaxter 
2010), analyses of genetic diversity, and phylogenetics (Cruaud et al. 2014). Additionally, RADseq is advantageous 
over whole-genome sequencing because it produces many independent amplifications of the target regions, yield¬ 
ing a greater depth of coverage, which is important to distinguish sequence errors from sequence variants (Sims 
et al. 2014), and because it is relatively inexpensive, allowing more samples to be sequenced without requiring 
prior information on the genomes of interest (Andrews et al. 2016). Although generally not considered feasible for 
higher-level, broad scale phylogenetics, RADseq has been successfully used to analyze family-level relationships 
(Longo & Bernardi 2015). 

Objectives. In this study, our objectives were, first, to infer the phylogenetic relationships of a large number of 
species of liparids, with multiple individuals when possible, using readily available and easily obtained sequence 
data from the mitochondrial COI gene. New sequence data were obtained from species collected and identified 
from the eastern North Pacific and augmented with data from species in the BOLD and GenBank public databases 
to provide insight into global interrelationships. Second, we attempted to infer a more robust phylogeny for a select 
group of eastern North Pacific species for which we could obtain an appropriate quality of data by using the more 
intensive techniques of RADseq analysis to survey hundreds of genome-wide restriction-site associated loci. 


Materials and methods 

Sample collection. Specimens were collected primarily during survey operations conducted in the east¬ 
ern North Pacific from Alaska to California by the Resource Assessment and Conservation Engineering 
(RACE) Division of the U.S. National Marine Fisheries Service, Alaska Fisheries Science Center (AFSC). Ad¬ 
ditional specimens from the Salish Sea were obtained from the field collections conducted by the Washing¬ 
ton Department of Fish and Wildlife and University of Washington. Tissue samples were taken from fresh 
whole specimens at sea or from frozen specimens in the laboratory and preserved in 95% ethanol for storage 
at -80°C. Most whole specimens serving as vouchers were fixed in 10% formalin and stored in 70% ethanol; 
some whole specimens were fixed and stored in 95% ethanol. Right-side pectoral girdles, from specimens fixed 
in 10% formalin and preserved in 70% ethanol, were dissected, cleared, and counter stained following Pot- 
thoff (1984). Institutional codes follow Sabaj (2016). All specimens and tissues were archived at the Burke Mu¬ 
seum of Natural Flistory and Culture, University of Washington Fish Collection (UW; Appendix Tables 1-3). 
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Careproctus notosaikaiensis 
Careproctus longipectoralis 
Careproctus continentalis 
Careproctus bowersianus 
Careproctus attenuatus (C. longifilis) 
Elassodiscus tremebundus 
Elassodiscus caudatus 2 
Elassodiscus caudatus 1 
Paraliparis rosaceus 
Paraliparis sp. 

Paraliparis bathybius 

Crystallichthys cyclospilus 

Crystallia matsushimae 

Careproctus georgianus (Paraliparis sp.) 

Paraliparis mawsoni 

Paraliparis neelovi 

Paraliparis sp. 

Paraliparis sp. (P. camilarus) 

Paraliparis paucidens 

Paraliparis aff. longipectoralis (P. macropterus) 
Paraliparis antarcticus 
Paraliparis sp. (P. antarcticus) 

Paraliparis charcoti 
Paraliparis leobergi 
Paraliparis sp. 

Paraliparis sp. 

Paraliparis valentinae 
Edentoliparis terraenovae 
Liparidae sp. (Paraliparis pectoralis) 

Paraliparis pectoralis 
Liparidae sp. (Paraliparis mento) 

Rhinoliparis barbulifer (R. attenuatus) 
Rhinoliparis sp. (R. attenuatus) 

Rhinoliparis attenuatus 

Paraliparis melanobranchius (P. megalopus) 

Paraliparis thalassobathyalis 

Paraliparis cephalus 

Lipariscus nanus 

Careproctus sp. (C. cypselurus) 

Careproctus furcellus (C. cypselurus) 
Careproctus cypselurus 
Careproctus sp. cf. melanurus 
Careproctus melanurus 
Careproctus sp. (C. furcellus) 

Careproctus canus (Careproctus n. sp.) 
Paraliparis sp. (Careproctus sp.) 

Careproctus sp. (C. colletti) 

Careproctus colletti 
Careproctus simus 
Careproctus ostentum 
Careproctus gilberti 
Careproctus rastnnus (C. scottae) 

Careproctus scottae 
Careproctus reinhardti 
Careproctus phasma 
Paraliparis dactylosus 
Careproctus comus 
Careproctus faunus 
Acantholiparis opercularis 
Liparidae sp. 1 
Liparis pulchellus 
Liparis tessellatus 
Liparis rutteri 
Liparis sp. (L. dennyi) 

Liparis dennyi 

Liparis mucosus (L. dennyi) 

Liparis ochotensis 
Liparis gibbus 
Liparis bathyarcticus 
Liparis tanakae 
Liparis chefuensis 
Liparis agassizii 
Liparis liparis 
Liparis montagui 
Liparis tunicatus 
Liparis inquilinus 
Liparis fabricii 
Liparis sp. 

Liparis atlanticus 
Liparis callyodon 
Polypera greeni 
Liparis florae 
Liparis fucensis 
Nectoliparis pelagicus 
Cyclopterus lumpus 
Eumicrotremus spinosus 
Eumicrotremus derjugini 
Aptocyclus ventricosus 


FIGURE 7. Majority-rule (50%) consensus phylogenetic tree of Shen et al. (2017, after fig. S6), derived from a Bayesian 
inference of a 440 bp alignment of cytochrome c oxidase subunit 1 gene (COI) sequences for 84 samples of 83 liparid species. 
Bayesian posterior probabilities are above branches. Tree is rooted with species of the Cyclopteridae. Corrected identifications 
based on our study are in parentheses. 
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Cytochrome c oxidase subunit I sequences. We attempted to collect sequence data for species never before 
analyzed as well as for species used in previous molecular phylogenetic or barcoding analyses. With the exception 
of Psednos, we included representatives of all valid genera and all species in these genera used in previous broad- 
scale analyses (i.e., Knudsen et al. 2007; Rock et al. 2008; Steinke et al. 2009; Duhamel et al. 2010; Gardner et al. 
2016; Shen etal. 2017). We used a total of 268 sequences of at least 122 species among 18 genera plus two outgroup 
taxa for the COI phylogeny (Appendix Table 1). In deciding which data to use, we evaluated over 773 sequences 
derived from 247 tissues of new specimens collected primarily from across the North Pacific and 526 sequences 
from the public online databases BOLD and GenBank (Appendix Tables 1 and 2). Only sequences among species 
that were complete and not identical were included in the dataset used in the phylogenetic analysis. Other sequences 
surveyed and identical to those used in the analysis are listed in the Appendix Table 2. Additionally, for clarity in fig¬ 
ures, monophyletic species clades were collapsed, except when notable sequence divergences were evident—these 
are discussed in the text. Illustrated branch tips are thus often supported also with specimens representing multiple 
haplotypes. Total numbers of sequences are in parentheses following taxa on the tree and each sequence is listed in 
Appendix Tables 1 and 2. 

We used sequences from public databases with special caution because misidentifications among snailfishes are 
common within them. These data were used only when we were able to confirm identifications by verifying sources 
from vouchered specimens or their photographs, by comparison with our new data, or by accepting published 
identifications. A complete list of material with corrected identifications is provided in Appendix Tables 1 and 2. 
Our corrections to identifications in historical trees (Figs. 1-7) are presented in parentheses following the names as 
originally published. Institutional abbreviations are those provided by Sabaj (2016), unless otherwise indicated. 

To obtain new sequence data, we extracted DNA using a QIAGEN DNeasy blood and tissue kit (Valencia, CA, 
USA) in single column format. We used primers HC02198 and LCO1490 to amplify a 710 bp segment of the COI 
gene via the polymerase chain reaction (PCR) in 25 pi volume, as described in (Folmer et al. 1994). PCR products 
were sequenced in both forward and reverse directions with Sanger sequencing using the PCR primers at the High 
Throughput Genomics Unit at the University of Washington (Seattle, WA, USA), with the exception of Progna- 
tholiparis ptychomandibularis (UW 116036), Careproctus Candidas (UW 156091-1, UW 156091-2), C. staufferi 
(UW 155711, UW 119196-1), Paraliparis penicillns (UW 119192), and P. ulochir (UW 119776-1), which were 
sequenced at MCLAB (San Francisco, CA, USA). Contigs were assembled, checked manually with their chromato¬ 
grams using Sequencher (2011, version 5.0, Gene Codes Corporation, Ann Arbor, MI), and aligned using BioEdit 
version 7.2.6 (Hall 1999). Forward and reverse sequences were obtained from all 247 new samples. A fragment of 
490 bp was used for analysis after trimming ends with low sequencing quality. New sequences were deposited in 
BOLD or GenBank. BOLD sequence identification numbers and GenBank accession numbers are given in Appen¬ 
dix Tables 1 and 2. 

Restriction-site Associated DNA sequences. We used tissues from 44 individuals of 28 species among 10 
genera of liparids, plus four individuals of the cyclopterid Eumicrotremus orbis, to generate the phylogenetic tree 
from RADseq data (Appendix Table 3). All sequences were generated from samples taken from whole specimens 
from the eastern North Pacific archived at the University of Washington (Appendix Table 3). DNA was extracted 
using the same methodology as for COI sequencing. DNA was quantified using the Quant-iT PicoGreen dsDNA 
Assay Kit (Invitrogen Inc., Carlsbad, CA), and GEN5 software (BioTek, Winooski, VT). DNA extraction produced 
low yields, eliminating some species (e.g., Liparis fucensis), which we suspect was due to the presence of excessive 
mucopolysaccharides in snailfishes (Jeong 2008). We designated species groups (Appendix Table 3) in the DNA 
library. Some species groups were represented by four individuals; many by only one. DNA of the selected samples 
was concentrated using a vacuum concentrator to 6.25 ng/pl in a 20 pi volume. Preparation of pooled RAD-tagged 
fragments (Sbjl restriction enzyme digestion, adapter ligation, shearing, and PCR) was performed following the 
methods of Baird et al. (2008) and Hohenlohe et al. (2011). Following DNA digestion with Sbfl, PI adapters were 
ligated and a unique 6-nucleotide barcode was assigned to each individual pooled into a single library. DNA was 
randomly sheared to a target size of 300-500 bp using a Qsonica Q800R DNA sonicator (Qsonica, LLC, Newtown, 
CT). Size selection was performed using Agencourt AMPure XP magnetic beads (Agencourt Bioscience Corp., 
Beverly, MA). A second set of adapters (P2) was then added using blunt-end ligation. The base-pair target length of 
libraries (200-400 bp), concentration, and quality were assessed with a 2100 Bioanalyzer and DNA 1000 kit (Agi¬ 
lent Technologies, Santa Clara, CA). Single-read next-generation sequencing was performed as a single library in 
one direction at 150 bp target length (SRI 50) on a HiSeq 4000 at the University of Oregon Genomics Core Facility 
(https://gc3f.uoregon.edu/). 
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The Stacks software package v. 1.44 (Catchen et al. 2011, 2013) was used to filter raw genetic data to distin¬ 
guish sequencing errors from nucleotide polymorphisms within samples and to identify orthologous loci across 
taxa. The restriction site overhang (6 bp) and unique barcode identifier (6 bp) were removed from raw reads, which 
were then trimmed to 94 bp (due to low quality scores on the 3’ end of reads) and quality filtered using the pro¬ 
gram process radtags. To avoid genotyping errors associated with variable coverage, individuals with fewer than 
100,000 filtered reads were dropped from downstream analyses while the maximum number of filtered reads per in¬ 
dividual was capped at 2,000,000. The program denovo_map.pl pipeline in Stacks was used to call loci in individu¬ 
als and then build a catalog of loci from across all individuals. The minimum stack depth (-m) was five, the number 
of mismatches allowed between loci within an individual (—M) was two, and the number of mismatches allowed 
between loci among individuals when building the catalog (—n) was 10. Datasets with various amounts of missing 
data were generated by running multiple iterations of the Stacks program populations by altering filter parameters 
for the minimum number of species (—p) and the percentage of individuals within each species (—r) required to re¬ 
tain the locus. Complete RADseq loci were exported for each individual and concatenated into supermatrices using 
FASconCAT-G (Kiick & Longo 2014) with the most frequent option (-c) used for generating consensus sequences. 
The RADseq dataset of quality filtered reads was submitted to the National Center for Biotechnology Information 
(NCBI) sequence read archive (SRA) at https://www.ncbi.nlm.nih.gov/sra under accession number SUB4307188; 
SRA sample numbers are provided in Appendix Table 3. 

Zoogeographic sampling. Our COI dataset included a wide geographic range of species but was dominated by 
samples from the eastern North Pacific (Appendix Tables 1 and 2), while our RADseq dataset (Appendix Table 3) 
was made up entirely of specimens from the eastern North Pacific. Samples with COI data from the North Pacific 
were primarily from verified vouchers, although some were from online databases and not verified, as were most 
from the North Atlantic. All sequence data from the Southern Hemisphere, including species of Careproctus and 
Paraliparis (with Edentoliparis ), were taken from online databases with accepted identifications from the literature, 
limiting a more detailed interpretation of our results. Of the more than 127 species with COI data, 53% of the species 
included in the analysis were from the eastern North Pacific and Bering Sea, 24% from the Southern Hemisphere, 
17% from the western North Pacific, and 10% from the North Atlantic. While about half the known species of the 
Northern Hemisphere were sampled, only about 12% of Southern Hemisphere species were available. Our North 
Atlantic material was similarly limited, but reflective of the lower diversity of the region, with 13 species represent¬ 
ing four of the six genera in the region— Liparis, Careproctus, Paraliparis, and Rhodichthys —and about half the 
species. The 21 western North Pacific species sampled represent nearly half of the described liparids from the region 
(Nakabo & Kai 2013). 

Phylogenetic inference. Phylogenetic analyses of both the COI and RADseq datasets were perfonned under 
both maximum likelihood and Bayesian optimality criteria. The maximum likelihood analyses were conducted 
in R version 3.4.1 (R Core Team 2017) using the packages APE 4.1 (Paradis et al. 2004) and phangorn 2.2.0 
(Schliep 2011); the Bayesian analysis was conducted with MrBayes v3.2 (Huelsenbeck & Ronquist 2001; Ronquist 
& Huelsenbeck 2003). The modelTest function in the phangorn package was used to determine the most appropriate 
nucleotide substitution model for both data sets. In both cases, the model with both the greatest log likelihood score 
and the lowest AIC score was a generalized time reversible (GTR) model with gamma-distributed rate variation 
across sites and invariant sites (GTR + G + I). 

For the COI data set, we included at least 122 liparid species plus 10 individuals identified only to genus in the 
analysis, as well as two cyclopterid outgroup species ( Eumicrotremus orbis and Cyclopterus lumpus) to root the 
trees (Appendix Table 1). The final data matrix included 270 sequences, 490 bp in length. For the maximum likeli¬ 
hood analysis, parameters of the GTR + G +1 model were optimized using the optim.pml function of the phangorn 
package with nearest neighbor interchange (NNI) rearrangement, and 1,000 bootstrap replicates of the data set 
were created using the bootstrap.pml function of phangorn. Bootstrap support was mapped to the optimized tree 
using the plotBS function of phangorn, and the tree was rooted using the root function of APE. For the Bayesian 
analyses, posterior probability distributions were generated by running four Markov chains, under the default of 
three heated chains and one cold. After 10 million generations, the average standard deviation of split frequencies 
was well below 0.01, potential scale reduction factor (PSRF) scores were either 1.000 or 1.001 for all parameters, 
and minimum effective sample size (ESS) values were all well over 200, indicating that the runs had converged. 
Sampling frequency was 1,000 generations. The initial 2,500 samples were discarded as bum-in, and the remaining 
7,500 samples were used to estimate tree topology and posterior probabilities. 
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For the RADseq data set, we included 25 liparid species in the analysis, as well as two specimens of the cy- 
clopterid outgroup species Eumicrotremus orbis (Appendix Table 3). The final data matrix included 40 sequences 
33,370 bp in length from 355 RADseq markers, each 94 bp after trimming and concatenation. Procedures for both 
the maximum likelihood analysis and the Bayesian analysis were the same as for the COI data set, except that sat¬ 
isfactory convergence diagnostics for the MrBayes run (as noted above) were achieved in 4 million generations. 
Thus, the initial 1,000 samples were discarded as bum-in, and the remaining 3,000 samples were used to estimate 
tree topology and posterior probabilities. 


Results 

COI sequences. High levels of diversity were found among individuals with confirmed identifications of Careproc- 
tus melanurus, C. phasma, E. caudatus, and P. pectoralis. In each of these “species” at least two distinct monophy- 
letic taxa are evident. Some of this diversity had been recognized in the past (Steinke et al. 2009; Orr et al. 2015; 
Gardner et al. 2016), while other examples were newly discovered here after being compared with correctly identi¬ 
fied specimens. Our reidentifications based on vouchered specimens, either by physical examination of specimens 
or of photos, are listed in Appendix Tables 1 and 2. 

With few exceptions, all species were reciprocally monophyletic (Figs. 8-11). While we reidentified many 
records from online databases, exceptions were likely because of simple misidentifications of specimens in on-line 
databases, or previously unidentified cryptic variation, for which we could not verify identifications. Previously 
unidentified cryptic variation may also be responsible. The exceptions to monophyly occur among five groups of 
species: 

(1) Within the clade of Liparis agassizii, L. tanakae, and L. chefuensis (Fig. 9), we verified the identification of one 
specimen of L. tanakae (FAKU 135637), which differed from all others in this group by at least two base pairs. 
While identified in BOLD as L. agassizii, L. tanakae, or L. chefuensis, all other sequences were identical, except 
for that of a specimen of L. tanakae (ABJF123), which alone differed by a single base pair. 

(2) All specimens of L. fabricii were recovered with L. bristolensis in a polytomy basal to other Arctic and North 
Atlantic species. 

(3) Some specimens of Careproctus reinhardti, C. lerikimae, and C. phasma were recovered in a polytomy (Fig. 
10). Among several individuals identified as these species, sequences are identical, although most specimens of 
each species differ from all other species. 

(4) One individual of Paraliparis “bathybius” collected from the Southern Ocean (ANTFI351-06) was nested in a 
clade phylogenetically distant (Fig. 11) from 11 individuals of P. bathybius collected from the North Atlantic. 
Although listed with a pending species identification of P. thalassobathyalis (BOLD, 5 Jan 2018), it was distant 
also from the P. thalassobathyalis used in our analysis. 

(5) Three specimens of Paraliparis neelovi identified by Duhamel et al. (2010; MNHN 2008-2630, -2631, -2632) 
are apparently not the same species as specimens identified by Stein (2012; NMNZ R42667 and R42304). 
These specimens differ by 6.5-6.7 % sequence divergence (Fig. 11) as well as in the number of proximal pec¬ 
toral radials (Duhamel et al. 2010; Stein 2012). 

Among species, uncorrected sequence divergences ranged from 0 to 21 %, with an overall mean of 11.5%. 
With the exceptions noted above, 0% divergence was found among species of Paraliparis recently described from 
the Ross Sea (Stein 2012): among P. nullansa and P. plicatus, P. ekaporus, P. epacrognathus, P. posteroporus, and 
P. mentikoilon (Fig. 11). No divergence was also found between P. camilarus, P. longicaecus, and the previously 
described P. stehmanni Andriashev 1986; between P. andriashevi Stein and Thomplcins 1989 and P. valentinae 
Andriashev and Neyelov 1984; and between Paraliparis sp. (FOAG559) and P. neelovi Andriashev 1982 (MNHN 
2008-2631) (Fig. 11). The lack of sequence divergence among these species, despite described morphological dif¬ 
ferences (Stein 2012), indicates that the COI gene may not be not useful to distinguish among them (Smith et al. 
2012). Eliminating these 0% divergences, the lowest divergence was 0.2%, which was characteristic of many with- 
in-species comparisons throughout the Liparidae, as well as among other species of the Ross Sea. In addition to the 
cases listed above, intraspecific divergences above 1% were from largest to smallest within Liparis bathvarcticus, 
Paraliparis pectoralis, Liparis dennvi, Lipariscus nanus, and P. cephalus. 
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Within genera, uncorrected sequence divergences ranged from 0.2 to 9.1 %, with a mean of 4.0% (Table 1). In 
genera with multiple species represented, divergences ranged from 0.8% in Notoliparis to 9.1% in Careproctus. In 
genera represented by single species with multiple specimens, divergences ranged from 0.2% in Ciystallichthys and 
Pseudoliparis to 0.8% in Lipariscus. In the three major genera, sequence divergence was highest within Careproc¬ 
tus (9.1 %), followed by Liparis (8.5%) and Paraliparis (6.9%). Eliminating species of Paraliparis placed in distant 
clades (i.e., P. dactvlosus, P. grandis, and P. penicillus ), the sequence divergence within Paraliparis was reduced 
to 6.6%. 

Among genera, uncorrected sequence divergences ranged from 0.5 to 17.3 %, with a mean of 10.5% (Table 1). 
The smallest divergence was between Pseudoliparis and Notoliparis at 0.5%, followed by Squaloliparis and Crys¬ 
tal lie hthys at 4.8%, and Rhodichthys and Paraliparis at 5.7%. The largest divergence was between Nectoliparis and 
Elassodiscus, followed by that between Nectoliparis and Ciystallichthys, Lopholiparis, Lipariscus, and Prognatho- 
liparis, all over 15% divergence. 

RAD sequences. The number of filtered reads per individual ranged from 2,597 to 38,478,039, although 
for downstream analyses, filtered reads were capped at 2,000,000 and individuals with fewer than 100,000 were 
dropped. Of 47 individuals originally sampled, seven failed to yield the minimum number of quality filtered reads 
(< 100,000) and were removed from further analyses: Careproctus bowersianus (UW 119291), C. candidus (UW 
156091-02), C. staufferi (UW 155711), C. melanurus (UW 150589), Ciystallichthys cvclospilus (UW 119187), 
Paraliparis dactvlosus (UW 152008), and Eumicrotremus orbis (UW 151223). This left 30 species groups for 40 in¬ 
dividuals (Appendix Table 3). Maximum likelihood trees were created using —p 26, —p 27, and —p 28 (the minimum 
number of species groups in which a locus must be present) and r = 0.5 (the minimum percentage of individuals in 
a species group required to retain the locus) with the denovojnap software subroutine in Stacks. The tree with the 
highest likelihood was —p 28 and —r 0.5, which was based on 355 phylogenetically informative RADseq markers 
with length of 94 bp each (for a total of 33,370 bp). 

Uncorrected mean sequence divergences among species ranged from < 0.01 to 3.1 %, with a mean of 1.3%. 
Several species exhibited > 0.01% divergences between individuals: A.jordani (0.19%), A. ungak (0.03%), C. sp. 
cf. melanurus (0.01%), C. phasma (0.02%), C. scottae (0.01%), and Elassodiscus sp. cf. caudatus (0.03%). 

Within genera for which multiple known species were included, uncorrected sequence divergences ranged from 
0.1 to 2.0 %, with a mean of 0.8% (Table 2). Sequence divergence was highest within Liparis (2.0%), despite its 
resolution as a monophyletic genus in this analysis (see below), followed by Paraliparis (1.4%) and Careproctus 
(0.44%), both of which were recovered as paraphyletic. Allocareproctus exhibited the lowest mean sequence diver¬ 
gence of 0.14%. 

Uncorrected sequence divergences among genera ranged from 0.8 to 2.9 %, with a mean of 1.7% (Table 2). The 
smallest divergence was between Lopholiparis and Careproctus, followed by Allocareproctus and Careproctus. The 
largest was between Liparis and Nectoliparis, followed by that between both genera and Lipariscus. 

TABLE 2. Summary of mean uncorrected pairwise sequence divergence estimates (in %), within and among genera, for 
the RADseq data for taxa included in this study (Appendix Table 3). Number of unique sequences in parentheses; NA = 
within-genus distance not calculated for genera for which only one specimen was sequenced. 
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Allocareproctus (4) 

0.1 











2 

Careproctus (15) 

0.4 

0.8 










3 

Ciystallichthys (2) 

0 

1.2 

1.3 









4 

Elassodiscus (4) 

<0.1 

1.0 

1.1 

1.2 








5 

Liparis (4) 

2.0 

2.3 

2.5 

2.4 

2.4 







6 

Lopholiparis (1) 

NA 

1.0 

0.8 

1.5 

1.3 

2.5 






7 

Paraliparis (3) 

1.4 

1.1 

1.1 

1.3 

1.2 

2.5 

1.3 





8 

Lipariscus (1) 

NA 

1.5 

1.7 

1.5 

1.6 

2.6 

1.8 

1.4 




9 

Nectoliparis (1) 

NA 

2.2 

2.3 

2.2 

2.2 

2.9 

2.5 

2.4 

2.6 



10 

Prognatholiparis (1) 

NA 

0.9 

1.1 

1.1 

0.9 

2.3 

1.2 

1.2 

1.5 

2.1 


11 

Rhinoliparis (1) 

NA 

1.4 

1.6 

1.3 

1.5 

2.6 

1.7 

1.4 

1.3 

2.5 

1.4 
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COI phylogeny. The 50% majority rule consensus tree resulting from the Bayesian analysis (Figs. 8—11) was 
characterized by a high level of support at shallow nodes but significant polytomies, especially among non -Liparis 
species, at deeper intermediate nodes. Nodes at the basal portions of the tree were well supported, with posterior 
probabilities > 0.9. Although nearly all nodes with high Bayesian posterior probabilities also had high bootstrap val¬ 
ues calculated in the maximum likelihood analysis, the significant exception being the basal node at the divergence 
of L. fucensis from other liparids, many clades that were recovered in the consensus tree were not well supported 
by bootstrap values. 

In a conservative assessment of this tree, the three major liparid genera— Liparis, Careproctus, and Paralipar- 
is —are paraphyletic. Liparis is represented by two clades, Careproctus by at least four, and Paraliparis by three or 
more, if Edentoliparis and Rhodichthys are recognized as valid. All other genera with multiple species or individuals 
are monophyletic, including Allocareproctus, Elassodiscus and Rhinoliparis, as well as the monotypic Lipariscus 
and Nectoliparis. Our phylogeny infers many relationships that differ from previously published discussions of 
liparid relationships and of names applied to species groups. Table 3 lists all taxa used in our molecular analysis 
and the names we have applied. We did not attempt to place all species of the Liparidae within the classification, 
although we include some discussion about other species when warranted. 
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FIGURE 8. Phylogeny of the Liparidae. Majority rule (50%) consensus tree from the Bayesian inference of a 490 bp align¬ 
ment of 270 cytochrome c oxidase subunit one gene (COI) sequences. Nodal values represent Bayesian posterior probabilities 
and bootstrap values from the maximum likelihood analysis (above and below, respectively). Species names are followed by 
a catalog number or BOLD “Sequence ID” number when represented by a sequence from a single specimen in our dataset. N 
indicates number of sequences, when multiple sequences support a branch tip. Only unique sequences were subjected to the 
analyses (Appendix Table 1); other identical sequences surveyed are listed in Appendix Table 2. Clades Liparis , Aenigmoliparia, 
and Paraliparia are depicted in Figures 9, 10, and 11, respectively. 

The inferred phylogeny from the COI dataset included a clade of individuals of Nectoliparis pelagicus recov¬ 
ered as the sister of all other liparids (Fig. 8). Nectoliparis is followed by Liparis fucensis, sister of all other liparids 
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and with the consequence of rendering Liparis paraphyletic, a node with a posterior probability of 1.0 but with low 
bootstrap support. Liparis fucensis is followed by a clade of all remaining species of Liparis (Clade A). Among these 
species of Liparis are four major clades (Fig. 9). A clade of L. florae + L. greeni + L. callyodon + L. curilensis is 
sister of all other Liparis. A polytomy follows containing a clade of L. punctulatus + L. rutterv, a clade of L. muco- 
sus + L. tessellatus + L. pulchellus', L. cyclopus in a separate lineage; and a fourth clade containing all remaining 
Liparis. The three smaller clades of this polytomy all consist of North Pacific species. The large remaining clade of 
this polytomy includes L. ochotensis alone as the sister of two larger clades. One includes several Arctic and North 
Atlantic species {L. fabricii, L. liparis, L. tunicatus, L. inquilinus, L. bristolensis, L. atlanticus, andT. montagui)', 
the other is composed of a polytomy with a clade of western North Pacific species (L. tanakae, L. agassizii, and 
L. chefuensis)', a clade of L. bathyarcticus; and a clade of L. dennyi + L. gibbus. While most relationships are well 
supported, the positions of L. mucosus, L. montagui, L. greeni, L. dennyi, and L. fabricii are not, and monophyly of 
both L. dennyi and L. fabricii, although not rejected, is not strongly supported. 
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FIGURE 9. Phylogeny of the genus Liparis, excluding L. fucensis depicted in Figure 8, from the majority rule (50%) consen¬ 
sus tree from the Bayesian inference of a 490 bp alignment of 270 cytochrome c oxidase subunit one gene (COI) sequences. 
Nodal values represent Bayesian posterior probabilities and bootstrap values from the maximum likelihood analysis (above and 
below, respectively). Species names are followed by a catalog number or BOLD “Sequence ID” number when represented by 
a sequence from a single specimen in our dataset. N indicates number of sequences, when multiple sequences support a branch 
tip. Only unique sequences were subjected to the analyses (Appendix Table 1); other identical sequences surveyed are listed in 
Appendix Table 2. 
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Following the Liparis clade (Fig. 8), a clade of Careproctus crozetensis and the hadal Pseudoliparis + Notoli- 
paris (Clade B in Fig. 8) is followed by a clade of the eastern North Pacific C. ovigerus and Antarctic C. georgianus 
(Clade C in Fig. 8). Both clades are highly supported. 

A large clade (Clade D in Fig. 8) follows and includes the two major clades E and F. Clade E has both low poste¬ 
rior probability and low bootstrap support and is composed of a large polytomy of several small clades and the large 
and highly supported Clade G (Fig. 10). The polytomy comprises Acantholiparis, Allocareproctus, and Prognatho- 
liparis, C. longifllis, and Paraliparis penicillus; a clade of three Antarctic species of Careproctus with the eastern 
North Pacific Osteodiscus cascadiae; a clade of undescribed North Pacific Careproctus with C. canus; a clade of 
Elassodiscus with C. bowersianus, C. notosaikaiensis, and C. cyclocephalus; a clade of Careproctus candidus + 
C. comus + C.faunus + C. iacchus + C. Stauffer i\ and the large Clade G. Clade G is composed of Lopholiparis as 
the sister of three well-supported clades of Careproctus with a clade of two North Pacific species of Paraliparis ( P. 
dactvlosus and P. grandis) among them. Clade N contains the North Pacific C. colletti + C. cvpselurus + C.furcel- 
lus + C. melanurus plus a clade of two unidentified southern Careproctus with an undescribed (Orr, unpublished) 
northern Careproctus, and is the sister of a clade of the North Pacific C. simus and Arctic Atlantic C. micropus and 
C. kidoi. Clade O contains three lineages: the C. rastrinus species group (C. rastrinus + C. scottae + C. acanthodes 
+ C. trachysoma + C. pellucidus); the C. gilberti species group (C. ostentum + C. gilberti)', and the C. reinhardti 
species group (C. phasma + C. lerikimae + C. reinhardti). Clade P encompasses both P. dactvlosus and P. grandis, 
recovered in a polytomy with clades N and O. 

Clade F has high posterior probability and high bootstrap support and includes two well-supported clades. 
The first contains Ciystallichthys, Squaloliparis, and western Pacific species of Careproctus : C. macrodiscus, C. 
marginatus, and C. roseofuscus (Clade FI). The second clade (Clade Q), which is larger and more poorly supported, 
includes nearly all derived species that lack a pelvic disc (excepting only Acantholiparis opercularis, and three 
members of Paraliparis — P. penicillus, P. dactvlosus, and P. grandis). 

Within Clade Q, Paraliparis pectoralis is sister of all other Paraliparis, although the position is not well sup¬ 
ported (Fig. 11). Lipariscus, Rhinoliparis, and six clades of Paraliparis, including Rhodichthys and Edentoliparis 
embedded within, arise from a large polytomy. Rhinoliparis and Lipariscus are both monophyletic and well sup¬ 
ported. Other species of Paraliparis fonn five highly supported clades arising from the polytomy: P. wolffi + P. 
megalopus (reidentified from “ P. melanobranchus ”) + P. ulochir (Clade I); P. “neelovi ” + P. copei + P. rosaceus + 
P. “bathybius” (an Antarctic specimen likely misidentified) (Clade J); P. cephalus + P. mento (Clade K); Rhodich¬ 
thys + P. operculosus + P. charcoti + P. antarcticus + P. leobergi + P. valentinae + P. andriashevi + P. macropterus 
+ Paraliparis (= Edentoliparis) terraenovae (Clade L); and P. garmani + P. paucidens + P. mawsoni + P. bathybius 
plus recently described Paraliparis from the Ross Sea and unidentified Paraliparis of the eastern North Pacific 
(Clade M). A sixth lineage of Paraliparis is represented by P. thalassobathvalis. 

RAD phylogeny. Bayesian and maximum likelihood topologies of the RAD tree of the individual consensus 
sequences of eastern North Pacific species (Fig. 12) were identical. In contrast with the COI phylogeny (Figs. 8-11), 
all but three nodes of the RAD tree were well supported with Bayesian posterior probabilities of 1.0, and all but 
five nodes were supported with bootstrap values higher than 80. In this tree (Fig. 12), Nectoliparis was recovered as 
the sister of all other liparids. Liparis is monophyletic (although the analysis did not include L.fucensis, because of 
problematic extractions), followed by a clade that includes Ciystallichthys first, followed by a ladder arrangement 
of Rhinoliparis, Paraliparis ulochir, and Lipariscus. Prognatholiparis is the sister of a clade containing Elassodis¬ 
cus and Paraliparis penicillus, a monophyletic Allocareproctus, and a clade of Careproctus with its sister group 
Lopholiparis followed by a ladder arrangement of C. staufferi, C. candidus, C.faunus, and C. comus. Its sister clade 
comprises the Careproctus melanurus group, with the C. melanurus/C. sp. cf. melanurus clade sister of C. cypse- 
lurus and C. furcellus; Paraliparis dactvlosus; and the C. rastrinus group (C. scottae, C. gilberti, and C. phasma). 
The position of Lopholiparis has a high posterior probability (1.0) but low bootstrap support (76%). Although with 
posterior probability of less than one (0.97), the internal nodes of the Careproctus rastrinus clade are also well sup¬ 
ported. The position of Paraliparis dactylosus as sister of the C. rastrinus complex has a 1.0 posterior probability 
but among the lowest bootstrap support within the tree. Careproctus and Paraliparis are paraphyletic, with species 
of both genera found in three clades. Three nodes have a posterior probability < 1.0 and low bootstrap support: 
within E. caudatus, within C. melanurus, and between C. gilberti and the clade of two individuals of C. phasma. 
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Discussion 


General classification. Our phylogenetic results support the integrity of some elements of current classifications, 
although Liparis, Careproctus, and Paraliparis appear to be paraphyletic. While our results suggest a paraphyletic 
Liparis, the low support for the relationship leads us to reject recognizing multiple genera within Liparis without ad¬ 
ditional molecular or morphological evidence. Although as currently understood, Paraliparis is paraphyletic, a core 
component of species likely forms a monophyletic Paraliparis , with Edentoliparis and Rhodichthys in synonymy. 
However, three species of Paraliparis were recovered among species of Careproctus. In addition, other genera were 
uncertainly recovered within the major polytomy within Paraliparis. Careproctus is certainly paraphyletic but rec¬ 
ognizing multiple clades of the current Careproctus at the generic level is likely to lead to instability because of in¬ 
complete taxon sampling and the lack of support for many of the nodes that would lead to generic reassignments. 

Therefore, to encourage stability and to facilitate communication of information about clades we propose the 
classification detailed in Table 3 and illustrated in Figures 8M2. For the species we examined among the Liparidae, 
we used as appropriate the few published genus-group names and otherwise created new infonnal clade and species- 
group names for highly supported clades recovered in both our COI (Table 3) and RADseq analyses (Table 4). Other 
than allocating some species to available genus-group names, we propose no formal nomenclatural actions here. 

Paraphyly of Liparis. The paraphyly of Liparis with the inclusion of L. fucensis was a somewhat surprising 
result of our analysis. This arrangement has been reported in other analyses of COI data based on a single sample of 
L. fucensis from BOLD (Duhamel et al. 2010; Shen et al. 2017). We verified the identification of the voucher of that 
sample and included several additional individuals in our analysis, as well as many other species and individuals of 
Liparis (Appendix Tables 1 and 2). In our analysis, Liparis fucensis was monophyletic, with all specimens cluster¬ 
ing together, and was the sister group of all other liparids except Nectoliparis. However, bootstrap support for this 
paraphyletic relationship was weak and additional data should be incorporated before drawing strong taxonomic 
conclusions. No adult morphological characters have been identified that would suggest this primitive position out¬ 
side of other species of Liparis, and most recently Chernova (2008) placed it among other species in the subgenus 
Neoliparis, with L. bikunin, L. tarticus, and L. cyclopus, of which only L. cyclopus was included in our analysis. 
However, larvae of L. fucensis have a bubble-morph morphology (Marliave & Peden 1989), a derived character 
unknown in other species of Liparis. A similar, but different, larval morphology has been recorded in Nectoliparis 
pelagicus (M. S. Busby, pers. comm., 30 Jan 18) as well as in psychrolutid sculpins (Marliave & Peden 1989). 

Polypera Burke is a generic name originally erected for Liparis greeni (Burke 1912) and has also included Li¬ 
paris simushirae (Chernova et al. 2004; Parin etal. 2014). Recently, the genus has been considered valid (Chernova 
et al. 2004; Parin et al. 2014) or synonymous with Liparis (Kido 1988; Mecklenburg et al. 2002; Pietsch & Orr 
2015). In our analysis, Liparis greeni is placed within a well-supported clade that includes L. florae, L. curilensis, 
and L. callyodon, providing no support for the validity of the genus, at least when conceived as containing L. greeni 
alone. 

Subgenera of Liparis. Eight valid genus-group names have been proposed for Liparis (Chernova et al. 2004; 
plus Polypera), and Chernova (2008) recently reallocated species of Liparis among five of the subgenera. Of these, 
Liparis, Careliparis Garman, Lycocara Gill, and Lvoliparis Jordan and Evennann may accurately reflect aspects of 
the phylogeny of the genus; Neoliparis Steindachner is paraphyletic. 

Among the species allocated to the subgenus Liparis (Chernova 2008), with type species L. liparis, we included 
L. bristolensis, L. liparis, and L. tunicatus in our analysis. All were resolved in a polytomy with individuals of L. 
fabricii (Fig. 9), which had been placed in the subgenus Lycocara with L. koefoedia, the type species, and L. laptevi 
(Chernova 2008), for which we had no material. Also in the polytomy were species previously placed in Neoliparis 
(Chernova 2008; L. montagui, L. atlanticus, and L. inquilinus). Future resolution of this polytomy, by including 
presently omitted species or additional support from morphological or new molecular data, may recover a mono¬ 
phyletic subgenus Liparis that does not include Lycocara (L. fabricii). We recommend incorporating L. bristolensis, 
L. tunicatus, L. montagui, L. atlanticus, and L. inquilinus with L. liparis into the subgenus Liparis and retaining 
Lycocara for L. fabricii. 
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FIGURE 10. Phylogeny of the liparid clade Aenigmoliparia from the majority rule (50%) consensus tree from the Bayesian in¬ 
ference of a 490 bp alignment of 270 cytochrome c oxidase subunit one gene (COI) sequences. Nodal values represent Bayesian 
posterior probabilities and bootstrap values from the maximum likelihood analysis (above and below branches, respectively). 
Species names are followed by a catalog number or BOLD “Sequence ID” number when represented by a sequence from a 
single specimen in our dataset. N indicates number of sequences, when multiple sequences support a branch tip. Boldface spe¬ 
cies names indicate species placed in different positions in COI and RADseq trees. Only unique sequences were subjected to the 
analyses (Appendix Table 1); other identical sequences surveyed are listed in Appendix Table 2. 
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FIGURE 11. Phylogeny of the liparid clade Paraliparia from the majority rule (50%) consensus tree from the Bayesian infer¬ 
ence of a 490 bp alignment of 270 cytochrome c oxidase subunit one gene (COI) sequences. Nodal values represent Bayesian 
posterior probabilities and bootstrap values from the maximum likelihood analysis (above and below, respectively). Species 
names are followed by a catalog number or BOLD “Sequence ID” number when represented by a sequence from a single speci¬ 
men in our dataset. N indicates number of sequences, when multiple sequences support a branch tip. Only unique sequences 
were subjected to the analyses (Appendix Table 1); other identical sequences surveyed are listed in Appendix Table 2. 
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Careliparis, as defined by Chernova (2008), is not fully supported here. Rather, the seven species included in 
that subgenus formed three groups: L. ochotensis ; L. tessellatus (as sister of L. pulchellus ); and L. dennyi, L. gib- 
bus, and L. bathyarcticus plus a clade of L. agassizii (the type species), and L. tanakae. The latter five species form 
a well-supported clade with L. chefuensis (placed in Neoliparis by Chernova 2008; Fig. 9). We use the subgenus 
Careliparis in a more restrictive sense than Chernova and apply it only to this latter clade of six species. 
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FIGURE 12. Phylogeny of selected eastern North Pacific liparids inferred using genome-wide restriction-site associated DNA 
sequences (RADseq; -p 28, —r 0.5) with maximum likelihood and Bayesian methods. Majority rule (50%) consensus tree of in¬ 
dividual sequences. Nodal values represent Bayesian posterior probabilities and bootstrap values from the maximum likelihood 
analysis (above and below branches, respectively); double asterisks denote Bayesian posterior probabilities of 1 and bootstrap 
support of 100%. Species names are followed by the University of Washington Fish Collection catalog number for the specimen. 
Boldface species names indicate species placed in different positions in COI and RADseq trees. 
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Lyoliparis was erected by Jordan and Evermann (1896) as a monotypic subgenus for Liparispulchellus, a dis¬ 
tinctive elongate species with the caudal fin nearly entirely joined to the dorsal and anal fins. Liparis tessellatus is 
thought to have a convergent general morphology (Chernova 2008), but these data suggest rather that it is closely 
related to L. pulchellus despite some morphological differences (Chernova 2008). 

The subgenus Neoliparis as presently understood (Chernova 2008) is paraphyletic in our analysis, as species 
previously included in Neoliparis were placed among multiple clades. The type of the subgenus is L. mucosus, a 
species relatively unrelated to L.fucensis and other species assigned to Neoliparis (Chernova 2008); it alone of Neo¬ 
liparis is placed as the sister species of the Lyoliparis clade. In addition, Liparis fucensis was placed as the sister of 
all other liparids, rendering Liparis paraphyletic. Among other species, a clade of North Atlantic Liparis assigned to 
Neoliparis (Chernova 2008; L. atlanticus and L. montaguf) was distantly placed in a derived position among species 
of the subgenera Liparis and Lycocara. 

In our analysis, four species form a well-supported clade that is sister of all other Liparis, except L. fucensis : 
Liparis greeni, L. florae, L. callyodon, and L. curilensis (Fig. 9). We propose the subgenus Polypera, with type spe¬ 
cies L. greeni, for this clade. 

Paraphyly of Carepvoctus. Our analysis confirms that Careproctus is paraphyletic and will eventually require 
extensive taxonomic revision. However, because many intermediate nodes of our tree are polytomies with several 
clades of Careproctus mixed with clades of other genera, we choose to retain a paraphyletic Careproctus for con¬ 
venience and propose several informal names for clades, in some cases resurrecting older generic or subgeneric 
names. Also, four clades were composed of multiple genera including species allocated to Careproctus : the Pseu- 
doliparia (Fig. 8), with C. crozetensis the sister of Notoliparis and Pseudoliparis'. the Crystallichthyia (Fig. 8), with 
Ciystallichthys and Squaloliparis recovered within a clade of western Pacific Careproctus, and the Osteocareprocta 
(Fig. 10), with Osteodiscus recovered within a clade of southern Careproctus. For each of these clades, formal taxo¬ 
nomic treatments will be required when their relationships with other clades and genera are clarified with additional 
species and data in future analyses. 

Subgenera of Careproctus. Several genus-group names have been synonymized with Careproctus and are 
available. However, none precisely reflects our present understanding of liparid phylogenetics based on these re¬ 
sults. Two subgenera erected by Andriashev (2003) and further elaborated upon by Chernova (2005ab, 2014, 2015) 
have been recently diagnosed morphologically: Careproctus with 22 species and Careproctula with 46 species, each 
found in both Northern, nearly all in the Arctic and Atlantic oceans, and Southern hemispheres (Chernova 2005ab; 
Chernova 2014). Most species of Careproctus of the North Pacific have not previously been assigned to subgenera, 
which were often erected as genera or subgenera for single species at the time of their original description (e.g., Al¬ 
io chir Jordan & Evermann 1896 for C. melanurus, Prognurus Jordan & Gilbert 1898 for C. cypselurus, Temnocora 
Burke 1930 for C. candidus, Allinectes Jordan & Evermann 1898 for C. ectenes, and Caremitra Jordan & Evermann 
1896 for C. simus). 

We refer Clade O (including the members of the C. rastrinus species group of Orr et al. 2015; the C. gilberti 
species group; and the C. reinhardti species group) to the subgenus Careproctus (Figs. 10, 12). Chernova (2005ab, 
2014, 2015) assigned 14 northern species of the North Atlantic and Arctic to the subgenus Careproctus, with the 
northern type species C. reinhardti. Included in the subgenus are also seven southern species, following Andria- 
shev’s (2003) concept of the subgenus. In our analysis, all members of this clade are northern species; none of the 
seven southern species was available for this analysis. Chernova (2005ab) further argued for the division of northern 
members from southern members of the subgenus on the basis of the long lower pectoral-fin lobe (“long-fmned 
careprocts”; Chernova 2005a), among other characters, which also characterizes nearly all members of the C. rast¬ 
rinus and C. gilberti species groups. Among northern species of the subgenus Careproctus, Chernova (2005a) also 
informally recognized two or more groups, including the long-fmned “slit-eyed careprocts.” None of the species 
in these groups was available for our analysis. Chernova (2005b) also placed Careproctus furcellus (= Cvclopterus 
gelatinosus Pallas of Chernova 2005b) in the subgenus, although it is not closely related according to our analysis 
(Figs. 10, 12). 

Careproctus reinhardti likely represents a species complex as suggested by Chernova (2005b), who clarified the 
morphological characters of the type of C. reinhardti and the confusion promulgated by successive authors based 
on misidentifications. Orr et al. (2015) described C. lerikimae as a new species based on material from the Beau¬ 
fort Sea; sequences of some specimens are identical to that of other specimens identified as C. reinhardti. These 
specimens of C. reinhardti are likely misidentified, as paratypes of C. lerikimae had previously been identified as 
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C. reinhardti in the BOLD database (e.g., DSFIB 768-11, DSFAL600-09, DSFIB63-11; Orr et al. 2015; Appendix 
Tables 1—2). In addition, the illustration associated with C. reinhardti in Mecklenburg et al. (2016) is clearly not the 
species diagnosed by Chernova (2005b) and is more similar to several other species recently described by others 
(Chernova 2005ab, 2014; Orr et al. 2015). Nearly identical sequence data were extracted from other specimens that 
may represent additional species; and a specimen of C. phasma collected in the Gulf of Alaska (UW 154442) may 
represent a species different from C. phasma of the Bering Sea (Orr et al. 2015). 

We refer some members of Clade N (Clade Melanura) to the subgenus Allochir with type species C. melanurus 
(Jordan & Evermann 1896): C. melanurus, C. colletti, C. cypselurus, C.furcellus, an undescribed species originally 
labeled as “C. canus ” in BOLD from British Columbia, and an unidentified southern Careproctus. Two subgeneric 
names, plus subgenus Careproctus for C. furcellus, have previously been proposed for two species of this group: 
the senior name Allochir was erected for C. melanurus and the junior name Prognurus, for C. cypselurus. We refer 
other members of Melanura, previously considered as members of the subgenus Careproctula, to the subgenus 
Caremitra. 

The subgenus Careproctula Andriashev is thought to include 46 species that are each restricted to either the 
Northern or Southern hemisphere (Andriashev 2003; Chernova 2005b; Chernova 2014, 2015). The type species is 
C. fedorovi, a southern species that was not included in our material. Of the 46 species originally assigned to Care¬ 
proctula (Chernova 2014), our material included only five species. Three were among the 37 species of the south¬ 
ern group (Chernova et al. 2004; Chernova 2014): C. continentalis, C. longipectoralis, and C. georgianus. These 
southern species were resolved into two clades distant from the northern clade and from one another (Figs. 8, 10). 
The clade comprising the southern species C. continentalis, C. longipectoralis, and C. discoveiyae also includes the 
eastern Northern Pacific Osteodiscus cascadiae (C. cypselurus I and II of Steinke et al. 2009). To avoid creating a 
paraphyletic Careproctula, and rather than recognizing Osteodiscus cascadiae as a derived species of Careproctula 
or allocating all species of the clade to the genus Osteodiscus, we propose to use the informal name Osteocareprocta 
for this clade. Careproctus georgianus is a member of the Bathyphasmia clade with the eastern North Pacific C. 
ovigerus (Fig. 8), as discussed below. 

Nine other species originally assigned to Careproctula have been considered part of a northern Careproctula 
group (Chernova 2015). Of these, our material included Careproctus micropus and C. kidoi, which together formed 
a highly supported clade with C. simus within Clade Melanura (Fig. 10). For the clade containing C. simus, C. mi¬ 
cropus, and C. kidoi, we propose the resurrection of the subgenus Caremitra, with type species C. simus (Jordan & 
Evermann 1896). In their description of C. kidoi, Knudsen and Moller (2008) recognized that the species fit most 
closely with northern members of Careproctula based on several characters, but because of the presence of ribs 
and two hypural plates, their specimens did not fit the diagnosis, having more affinity with Careproctus. However, 
despite the common name “ribless Careproctus”, the diagnosis of Careproctula (Chernova 2005b) listed pleural ribs 
as “absent, rare rudimentary ribs present” and six of nine northern members of the subgenus have ribs (Chernova 
2015). The type species of the subgenus Caremitra, Careproctus simus, has rudimentary ribs (Orr, unpublished). 
Because our material included only five species of the 41 species previously allocated to Careproctula, the name 
Careproctula will likely apply to some monophyletic group of the remaining 36 species that includes C. fedorovi. 

We propose the name Bathyphasmia for the highly supported clade composed of Careproctus ovigerus of the 
North Pacific and C. georgianus of the Antarctic (Fig. 8). The genus Bathyphasma was erected by Gilbert (1896) for 
Careproctus ovigerus, primarily on the basis of its long, slender, sharp teeth and large pelvic disc. Balushkin (2012), 
in describing a new species of Volodichthys, argued for the resurrection of Bathyphasma based on a combination 
of general characters including canine teeth, a large gill slit, and an unnotched dorsal fin. However, Orr (2012) 
described two new species, both closely related to C. ovigerus on the basis of morphological characters alone and 
although very similar in other respects were distinguished from it by strongly trilobed teeth. The dorsal fin has a 
shallow notch in the smallest and largest specimens of C. ovigerus (Stein 1978; Orr 2012). Careproctus georgianus 
has a relatively derived pectoral girdle of four equally spaced radials lacking fenestrae (Andriashev 1998), unlike 
the plesiomorphic girdles of C. ovigerus, and its putative relatives C. lycopersicus and C. kamikawai (Orr 2012), 
as well as more vertebrae, simple teeth, and a small gill slit (Andriashev & Stein 1998). Because no known shared 
derived morphological characters diagnose the clade, we do not formally recognize the genus Bathyphasma. 

For Clade H (Fig. 8), composed of Careproctus macrodiscus, C. marginatus, C. roseofuscus, Crystallichthys 
cyclospilus, and Squaloliparis dentatus, we propose the name Crystallichthyia. Balushkin (1996) also recovered 
this clade in his analysis and noted the plesiomorphic condition of the unfused caudal fin retained by Careproctus 
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roseofuscus, C. marginatus, and S. dentatus, and characteristic also of Crystallichthys cyclospilus (Kido 1988), not¬ 
ing that further work is needed to identify morphological synapomorphies of the group. 

Within Clade E (Fig. 10), we recognize the subgenus Temnocora with type species C. candidus for the clade 
comprising Careproctus comus, C. faunus, C. iacchus, C. staufferi, and C. candidus. The genus Temnocora was 
originally erected for C. candidus alone because of its anterior dorsal-fin lobe and slit-like pupil (Burke 1930). 
While only C. candidus among these species has a distinctly notched dorsal fin, all except C. comus and C. faunus 
have similar slit-lilce pupils. Although recovered in the COI tree as the sister species of all other members of the 
Paracareprocta, Lopholiparis flerxi may be closely related, with its lobed dorsal fin and slit-like pupil, having been 
recovered as the sister of the Temnocora clade in the RADseq analysis, albeit with only moderate bootstrap support 
(Fig. 12). Tissues of other “slit-eyed” species of Chernova (2005a, 2014, 2015) were not available to us but these 
species are likely more closely related to members of the subgenus Careproctus. 

Paraphyly of Paraliparis. Our analysis confirms that Paraliparis is paraphyletic and will require some taxo¬ 
nomic revision. With the type species P. bathybius, a core derived component of species forms a monophyletic 
Paraliparis, and if extended to the clade including Edentoliparis and Rhodichthvs would place these genera in 
synonymy. Otherwise, Rhodichthys would be available for the clade including these genera. 

Within clade Paraliparia, Paraliparis terraenovae and Rhodichthys regina are both placed within a well-sup¬ 
ported clade of southern Paraliparis (Fig. 11: Clade L). Although P. terraenovae was allocated to the monotypic 
genus Edentoliparis by Andriashev (1990), a genus also later considered valid by him and others (Duhamel 1992; 
Andriashev 1998, 2003; Chernova et al. 2004), the genus was synonymized by Duhamel et al. (2010) on the basis 
of COI data, a conclusion supported as well by our analysis. The validity of Rhodichthys had been questioned by 
Andriashev (1986) and Kido (1988) but was reaffirmed with the description of a new species in the genus (An¬ 
driashev & Chernova 2011) and not examined by Duhamel et al. (2010). Based on our results, Rhodichthys should 
be synonymized with Paraliparis. Alternatively, if Rhodichthys is retained as a distinct genus, a new name will be 
required to accommodate species of the P. antarcticus species group. 

We are unable to firmly place Rhinoliparis and Lipariscus precisely among clades of Paraliparis because of a 
lack of resolution in the COI tree and lack of additional species of Paraliparis in the RADseq tree. Both genera were 
synonymized with Paraliparis by Kido (1988) because they were recovered among species of Paraliparis in his 
cladogram of 60 species (Fig. 1), similar to our results. However, in our COI analysis, these genera are uncertainly 
placed in the large polytomy with most other clades of Paraliparis (Fig. 11). Thus, our evidence neither supports 
nor explicitly rejects their synonymization with Paraliparis as proposed by Kido (1988) and, conservatively, we 
retain their generic status. 

The confirmation in this molecular analysis of the position of the Paraliparis dactvlosus species group (P. 
dactylosus and P. grandis ) among clades of Careproctus and other species with a complete pelvic disc (Fig. 10) is 
counterintuitive. However, earlier workers considered the species atypical of the genus (Stein 1978; Andriashev 
1986), and Stein (1978) stated that P. dactylosus was “not closely related” to species of Paraliparis known at the 
time. Paraliparis dactylosus may be unique among species of Paraliparis in having strongly trilobed teeth, although 
P. challengeri also has some teeth that have “well-developed lateral shoulders” (Chernova 2003, unnumbered fig¬ 
ure). In addition to its trilobed teeth, it has high counts of pyloric caeca, like many species of Careproctus. Most 
specimens of P. dactylosus we examined also clearly have a pair of pleural ribs, unknown in other Paraliparis 
except P. grandis. Paraliparis grandis itself was considered an “aberrant” Paraliparis by Andriashev (1986), and 
he eliminated several characters of the species from his generic description, including its unique pectoral fin and its 
high counts of pectoral-fm rays, pyloric caeca, and vertebrae. In addition, it has two pairs of pleural ribs, which are 
absent in nearly all other species of Paraliparis (Kido 1988; Stein et al. 2001). Among these characters, the pres¬ 
ence of pleural ribs, as well as high counts of pectoral-fm rays and pyloric caeca, are more typical for Careproctus 
than for Paraliparis. While many species of Careproctus bear the plesiomorphic condition of relatively strong ribs 
(Andriashev 1998; Chernova 2014, 2015), species of clade Melanura have one or two pairs of small, slender ribs 
(Orr, unpublished), unlike most species of the subgenus Careproctus, in which species of the C. rastrinus and C. 
reinhardti species groups have two or three strong, saber-like ribs (Chernova 2005b; Orr et al. 2015). Therefore, 
we suggest that Paraliparis dactylosus and P. grandis should both be considered derived species of Careproctus, 
closely related to the Melanura clade, and characterized by the independent loss of the pelvic disc. 

Paraliparis penicillus was recovered in the basal polytomy of clade Aenigmoliparia in the COI analysis (Fig. 
10) and sister of Elassodiscus in the RADseq analysis (Fig. 12). In both analyses it was distant from P. ulochir, the 
only other species of Paraliparis included in the RADseq data. 
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Subgenera of Paraliparis. As Stein noted in 1978, none of the subgenera of Paraliparis erected by earlier au¬ 
thors reflect the relationships as he understood them and with few exceptions (e.g., Chernova & Prut’ko 2011) later 
authors have declined to name groups of species. This uncertainty is evident in our analysis. Within Paraliparis, 
seven valid genus-group names have been formally proposed (Chernova et al. 2004). Of these we consider only the 
nominotypical subgenus Paraliparis and subgenus Amitrichthys to be appropriate for clades resolved in our study. 

The subgenus Paraliparis (Fig. 11: Clade M) includes species from nearly the entire geographic range of the 
Liparidae. The type species P bathybius (Collett 1879) and P. garmani are from the North Atlantic. The 15 other 
species comprising the clade include P. paucidens and an undescribed species from the eastern North Pacific, and 
several southern species, including those recently described from the Ross Sea (Stein 2012). 

Jordan & Evermann (1896) erected Amitrichthys for the type species P. cephalus. We propose to use this name 
for the clade composed here only of P. cephalus and P. mento (Fig. 11, Clade K). Paraliparis mento was also in¬ 
cluded by Jordan & Evermann (1896) in Amitrichthys, along with the distantly related P. rosaceus, P. copei, and P. 
dactylosus. 

Other names available for use as subgenera of Paraliparis include Amitra Goode 1881 (Monomitra Goode 1883 
is an unneeded replacement name; Eschmeyer et al. 2017) and Hilgendorfia Goode & Bean 1896. Amitra was a 
genus erected for Paraliparis liparinus, but the characters of P. liparinus are poorly known, based on a scant original 
description and an inadequate illustration of the holotype (Andriashev 1998; Chernova 2003). No other species have 
been proposed as closely related other than the almost equally rare P. challenged, and neither species was included 
in our analysis. The genus Hilgendorfia was erected for the type species P. membranaceus Gunther 1887 from off 
Chile. Although the original description and illustration were excellent, only the holotype is known and one other 
species from off Chile may be related to it (Stein 2005). Neither were available for our analysis. While Jordan & 
Evermann (1896) later allocated P. ulochir to Hilgendorfia, they provided no evidence to justify its closer relation¬ 
ship to P. membranaceus than to any other species. 

We propose the following informal species group names for other highly supported clades of Paraliparis. The 
Paraliparis copei species group (Chernova & Prut’ko 2011) includes/! copei and P. rosaceus, as well as specimens 
questionably identified as P. neelovi and P. “bathybius” from the Antarctic (Fig. 11, Clade J). The Paraliparis 
antarcticus species group includes P. antarcticus, P. macropterus, P. leobergi, P. charcoti, P. operculosus, Parali¬ 
paris (= Rhodichthys) regina, P. valentinae, P. andriashevi, and P. terraenovae (Fig. 11, Clade L). The Paraliparis 
dactylosus species group includes P. dactylosus and P. grandis (Fig. 10, Clade P) and appears to be distantly related 
to the other species of Paraliparis. The Paraliparis ulochir species group includes P. ulochir, P. megalopus, and P 
wolffi (Fig. 11, Clade I). 

Consensus between COI and RAD trees. Both COI and RADseq analyses supported similar topologies and 
conclusions about the phyletic integrity of major taxa. Allocareproctus was monophyletic. Careproctus and Parali¬ 
paris were paraphyletic, each split into at least three clades. Nectoliparis was the sister of all other liparids, followed 
by species of a monophyletic Liparis (except for L.fucensis, which was not available for the RADseq analysis). 

Significantly, both COI and RADseq analyses support the position of the Paraliparis dactylosus species group 
within and among otherwise highly supported clades of Careproctus. We eliminated the possibility of misidentifi- 
cation of these species by using multiple specimens that either shared the same COI haplotype or differed slightly 
but were placed within the same clade. In both trees, the relationship of the clade to either subgenus Allochir or 
subgenus Careproctus is unresolved—its position as sister of the subgenus Careproctus is marginally supported in 
the RADseq tree and it is in the polytomy With Allochir and subgenus Careproctus in the COI tree. 

Also, the fully resolved RADseq tree placed Crystallichthys as the sister of a clade comprising the disc-less 
Rhinoliparis, Lipariscus, and one member of Paraliparis. In the COI tree, Crystallichthys is placed within a highly 
supported clade of western Pacific Careproctus and Squaloliparis, which is sister of a large polytomy of clades of 
nearly all disc-less species we analyzed—a position that is highly supported. 

Incongruence between COI and RAD trees. The makeup of more derived clades and the relationships of 
several less diverse or monotypic taxa differed primarily because of the lack of resolution in the COI tree or the 
lack of comparative samples in the RADseq tree. Four clear conflicts are evident between the trees. Lopholiparis is 
sister of the Temnocora clade, composed of C. staufferi, C. candidus, C. comus, and C.faunus in the RADseq tree 
(Fig. 12); in the COI tree (Fig. 10), it is sister of the large, primarily northern, Clade Paracareprocta that includes the 
subgenera Allochir and Careproctus and the P. dactylosus species group. Its position is highly supported in the COI 
tree and is at least moderately supported in the RADseq tree. Within Temnocora, C. staufferi is sister of all other 
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members in the RADseq tree, while it is in a derived position as sister of C. candidus in the COI tree. In another ap¬ 
parent incongruity, Careproctus gilberti was placed between C. scottae and C. phasma in the RADseq tree (Fig. 12), 
but in a separate clade with C. ostentum in the COI tree (Fig. 10); RADseq data were not available for C. ostentum. 
The position of C. gilberti is strongly supported in both trees. Finally, Paraliparispenicillus is recovered in a clade 
with Elassodiscus caudatus in the RADseq tree (Fig. 12); but in the COI tree (Fig. 10), it is in the large polytomy of 
Clade Aenigmoliparia outside the Elassodiscus clade, within which P. penicillus is replaced with C. bowersianus, 
a species unavailable for RADseq. 

In addition to incongruences due to the lack of comparable samples, other differences in positions are apparent 
due to lack of resolution in the COI tree. In particular, Prognatholiparis is the sister of the remaining members of 
the Liparidae in the RADseq tree (Fig. 12), including an Elassodiscus/Paraliparis penicillus clade; in the COI tree 
(Fig. 10), it is in the large polytomy of Clade Aenigmoliparia with clades of Careproctus, Allocareproctus, the Elas¬ 
sodiscus clade (with C. bowersianus ), and P. penicillus. Similarly, Allocareproctus is in the same general position 
in both trees (Figs. 10, 12), although in the same large polytomy of the COI tree. 

Comparisons with previous results. Significant published phylogenetic hypotheses are those of Kido (1988; 
Fig. 1), Balushkin (1996; Fig. 2), Knudsen et al. (2007; Fig. 3), Steinke et al. (2009; Fig. 4), Duhamel et al. (2010; 
Fig. 5), Gardner et al. (2016; Fig. 6), and Shen et al. (2017; Fig. 7). In all analyses, all genera were monophyletic 
except for the three large, diverse genera Liparis, Careproctus, and Paraliparis. Liparis was considered monophy¬ 
letic until recent molecular analyses included L.fucensis ; the monophyly of Careproctus and Paraliparis has been 
suspect in most modem analyses. Our results further support the conclusion that these three genera are paraphyletic 
and require significant revision. 

The first modem phylogenetic analysis was conducted by Kido (1988; Fig. 1), who came to the surprising 
conclusion that the highly derived Nectoliparis pelagicus was the sister species of all other liparids. In all previous 
analyses, this pelagic snailfish—with its apomorphic characters of pelvic disc absent, gill slit entirely in front of the 
pectoral fin, and anus in a far anterior position—had been considered nested deeply within the family. Our analyses, 
as well as the results of other recent authors (Steinke et al. 2009; Gardner et al. 2016; Shen et al. 2017), strongly sup¬ 
port Kido’s (1988) conclusion, indicating that although highly derived, Nectoliparis or its ancestor diverged from 
the liparid lineage relatively early. In contrast, Balushkin (1996) argued for a more derived position for Nectoliparis, 
as traditionally accepted by previous authors (e.g., Burke 1930). 

Other than the position of Nectoliparis, Balushkin’s (1996; Fig. 2) phylogenetic hypothesis exhibits some broad 
similarities with the results presented here. His results, based on morphological data, include a basal clade of some 
species of Careproctus and three other genera, including Crystallichthys, in a polytomy comparable to our clade 
Crystallichthyia. In a second more derived clade, nine genera and other species of Careproctus are placed in a poly¬ 
tomy basal to a more derived clade comprising more species of Careproctus plus Elassodiscus, comparable to our 
clade Elassodisca. Finally, Balushkin (1996) recovered two derived clades of genera united by a highly reduced or 
absent pelvic disk: an unresolved component of eight genera, including the putatively monophyletic Paraliparis, 
and a clade of Nectoliparis and Lipariscus. Except for the derived position of Nectoliparis and Lipariscus, a mono¬ 
phyletic Paraliparis, and the inclusion of Acantholiparis, the clade is comparable to our Paraliparia. 

The first phylogenetic analysis incorporating molecular data was that of Knudsen et al. (2007; Fig. 3), who 
performed a combined analysis of mtDNA and morphological data for 24 species of the North Atlantic and North 
Pacific. In their final tree, Nectoliparis pelagicus was resolved as the sister of the rest of the Liparidae, followed by 
a monophyletic Liparis (again without L. fucensis), a monophyletic Careproctus, and a paraphyletic Paraliparis. 
Their limited taxon sampling included species of Careproctus of only our Paracareprocta clade and species of Para¬ 
liparis from two distinct clades: the subgenus Paraliparis and the P. copei species group. 

All subsequent analyses of liparid phytogeny (Rock et al. 2008; Steinke et al. 2009; Duhamel et al. 2010; Gard¬ 
ner et al. 2016; Shen et al. 2017) have been strictly molecular, using mitochondrial DNA sequence data, and in all 
cases both Careproctus and Paraliparis have been resolved as paraphyletic. Studies that included Liparis fucensis 
(Duhamel et al. 2010; Shen et al. 2017) have also resolved Liparis as paraphyletic, with L.fucensis as sister of all 
other Liparis. Our results concur with these findings, providing more evidence that L.fucensis may be inappropri¬ 
ately placed in Liparis. 

Paraliparis dactylosus has been placed by previous authors in positions different from that presented here. This 
apparent incongruence is the result of misidentifications of BOLD specimens (Appendix Tables 1 and 3). Steinke et 
al. (2009; Fig. 4) placed P. dactylosus with P melanobranchus (incorrectly labeled as “melanobranchius”; reidenti- 
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fied as P. megalopus) together in a clade with other species of Paraliparis and Lipariscus. Duhamel et al. (2010; Fig. 
5) recovered P. dactylosus with P. melanobranchus (also as “ melanobranchius ”) within a large polytomy comparable 
to our clade Paraliparia of disc-less species. Flowever, the sequence data for this “ dactylosus ” in both publications 
matches our new sequences of P. ulochir, and the photo published in BOLD is readily identifiable as P. ulochir. Our 
analysis recovered the clade of P. ulochir in our clade Paraliparia, comprising most species of Paraliparis, similar 
to the placement of “dactylosus” in both previous analyses. Shen et al. (2017; Fig. 7) used a correctly identified P. 
dactylosus (UW 119893) that was placed in a similar position to ours—among species of Careproctus. 


Morphology 

Among significant characters used for broad-scale phylogenetics in liparids are the absence, presence, and size of 
the pelvic disc; the detailed morphology of the pectoral girdle; and tooth patterns. Many other characters, such as 
cephalic pore patterns, meristics, and general morphometries, are critical for taxonomic analysis and are likely to 
prove useful for phylogenetic work, especially in more fine-scale analyses of the clades we identified (Stein 1978; 
Kido 1988; Knudsen et al. 2007). 

Pelvic disc. Perhaps the most significant character of liparid fishes is the highly modified pelvic fin, which is 
modified into a sucking disc or lost (Table 3). Our results suggest that the disc was independently lost at least three 
times: in the evolution of Nectoliparis, among members of the large disc-less clade of Paraliparis and other genera, 
and with the evolution of the Paraliparis dactylosus clade (Table 3). It was likely lost at other points as well, for 
instance in Acantholiparis and Paraliparis penicillus, both of which are placed uncertainly in our tree. When not 
lost, the disc evolved from the plesiomorphic condition of a large, well-developed disc similar to that found in the 
cyclopterids and the genus Liparis to typically becoming further reduced in more derived groups, as for example 
to the tiny, but complete, structure in Careproctus ostentum and C. gilberti. The disc was independently reduced to 
obsolescence in Elassodiscus. The unique skeletal pelvic disc of Osteodiscus (Stein 1978) is clearly a derived condi¬ 
tion from the large pelvic disc of related species. 

Pectoral girdle. The number and size of proximal radials, and notches, fenestrae, and the position of each have 
been considered important systematic characters of the pectoral girdle (e.g., Andriashev et al. 1977; Kido 1988; 
Andriashev 1998; Andriashev & Stein 1998; Stein et al. 2001). Distal radials, first identified in liparids by Orr & 
Maslenikov (2007), are under-reported but their presence and distribution may prove phylogenetically as well as 
taxonomically informative (Maslenikov et al. 2015). The plesiomorphic condition is a robust girdle with four proxi¬ 
mal radials (henceforth, “radials”), notched and relatively evenly spaced in the basipterygium, separated by three 
large fenestrae dorsally, a condition found in cyclopterid outgroups (Ueno 1970; Maslenikov et al. 2015; Fig. 13A, 
Table 3) and all species of Liparis (Andriashev 1986; Andriashev & Stein 1998; Fig. 13C). The more derived char¬ 
acter states show a reduction in all elements, ultimately to two round widely separated radials, either one dorsal and 
one ventral or both dorsal, without fenestrae (Andriashev 1990; Andriashev & Stein 1998; Table 3) or, uncommonly, 
enlarged evenly spaced radials lacking fenestrae (e.g .Elassodiscus and Acantholiparis, Fig. 131). Our results support 
the idea of independent and convergent reduction in pectoral girdle elements as first hypothesized by Andriashev 
(1998) for southern species. Although these ideas of the evolution of pectoral girdle morphology are supported in 
most respects by relationships based on molecular data, pectoral morphology is clearly not as tightly constrained as 
formerly proposed, and intraspecific variability among pectoral girdle characters is more common than previously 
expected (Knudsen & Moller 2008; Stein 2012, 2016). The plesiomorphic morphology is present basally in most of 
the major clades more derived than Liparis, including the deep-water clades Pseudoliparia, Bathyphasmia, Ambi- 
disca and its sister group Crystallichthyia, clade Aenigmoliparia and its major clade Paracareprocta. The discussion 
below pertains nearly exclusively to those species included in our analysis. A more complete picture of pectoral 
girdle evolution will become clearer with the incorporation of additional species into future phylogenetic analyses. 

The sister group of all other members of the Liparidae is Nectoliparis pelagicus, a highly apomorphic pelagic 
species. It possesses a pectoral girdle with three rounded radials in the upper half of the basipterygium while lacking 
the ventral fourth radial (Fig. 13B)—most similar to derived species of Paraliparis, as in P. ulochir (Fig. 130) and 
especially 12 Australian species described by Stein et al. (2001) for which we had no material. Most other species 
without a ventral radial have one or two dorsal radials (Kido 1988; Andriashev 1986, 1998; Stein et al. 2001). 

Among the hadal and bathyal Notoliparis, Pseudoliparis, and C. crozetensis, the pectoral girdle morphology of 
all except N. steward includes four round moderately sized radials without fenestrae (Duhamel et al. 2016; Stein 
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Careproctus georgianus Ant large R(l+1+1+1) simple Andriashev 1998 
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TABLE 3. (Continued) 

Lettered clade Taxa Distribution Pelvic disc size Pectoral radial formula Tooth shape, pattern Source 

Osteodiscus cascadiae eNPac large R(l+0+l) simple, bands RBCM 15941; Kido 1988; 

Stein 1978 
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Paraliparis antarcticus Ant absent R(2-3+l) simple, bands Andriashev 1986; Duhamel 

et al. 2010 

Clade M subgenus Paraliparis 

Paraliparis sp. Ant absent unknown unknown ANTFI326-06 

. continued on the next page 
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2016; Gerringer et al. 2017a). Notoliparis stewarti, recovered in a sister group relationship with P. swirei, has only 
three round radials (Stein 2016). Our results suggest that species of Notoliparis Andriashev 1975 may be derived 
members of Pseudoliparis Andriashev 1955, as suggested by Gerringer et al. (2017a). Other hadal species and gen¬ 
era will be needed to resolve this question conclusively. 

In the Bathyphasmia clade including C. ovigerus and C. georgianus, pectoral girdles are strikingly different. 
Careproctus ovigerus has a relatively plesiomorphic morphology with three deeply notched radials, with three large 
fenestrae, separated from the round ventral radial (Orr 2012), while C. georgianus has four round evenly spaced 
radials without fenestrae (Andriashev & Stein 1998), similar to those of the hadal and bathyal clade. 

Within the Ambidisca, among members of the clade Crystallichthyia, the sister group of the disc-less Paralipar- 
ia that also includes western Pacific species of Careproctus, the pectoral girdle morphology is of the plesiomorphic 
condition of three dorsal radials with one ventral radial. In all but C. macrodiscus, in which all radials are round and 
fenestrae are absent, the dorsal radials are notched and three fenestrae are present (Figs. 13D-G). 

Among clades of the Paraliparia, pectoral girdles are simplified and typically reduced. Fenestrae are rare or 
small when present. Exceptionally, among our included species, Paraliparis pectoralis (Fig. 13N) and P. stehmanni 
(Andriashev, 1986) have two small to moderately large fenestrae, and P. camilarus, Rhodichthvs melanocephalus, 
and Rhinoliparis attenuatus (Fig. 13L) have one small fenestra. Radials are typically round and notched only in 
those species with fenestrae, or uncommonly notched without fenestrae as in P ekaporus, P. epacrognathus, P. 
longicaecus, P. nullansa, and P. posteroporus (Stein 2012)—all members of the subgenus Paraliparis. The likely 
misidentified “P. neelovi” of Duhamel et al. (2010) placed in the P. copei species group has three dorsal radials, 
like P. copei (Andriashev 1986) but unlike P. neelovi of Stein (2012) which has only two dorsal radials, a more 
derived morphology congruent with this phylogeny. Among the four multiple-species clades for which pectoral-fin 
morphology is known, all members of the P. copei species group have three dorsal radials, while members of the 
P. ulochir species group all have a derived pattern with two dorsal radials (Fig. 130). Among members of the P. 
antarcticus species group, most have four radials evenly spaced and lacking fenestrae, except for the apomorphic P. 
terraenovae and Rhodichthvs (= Paraliparis) spp., which both have two dorsal and one ventral radials, and P. char- 
coti that uniquely among the species in our analysis has five radials (Matallanas 1999; although later described with 
only one or two dorsal radials and one ventral radial by Duhamel et al. 2010). Members of the subgenus Paraliparis 
have a range of two or three dorsal radials with one ventral radial, as well as some species having notched radials 
and fenestrae as noted above. In addition, Lipariscus nanus and P. thalassobathyalis, both with a highly reduced 
pattern of two dorsal radials and no ventral radial, and Rhinoliparis and P. cephalus (Fig. 13M), with the modified 
simple pattern of four large and equally spaced radials, are placed uncertainly. Further resolution of the Paraliparia 
polytomy with additional data will likely show that reduced girdles were independently derived from modifications 
of the plesiomorphic state of four equally spaced radials (Andriashev 1998). 

Among clades of the Aenigmoliparia within the large basal polytomy, the pectoral girdle morphology is un¬ 
known for several species. In some, a more plesiomorphic condition is evident. In Alio careproctus the girdle in all 
species is very similar to that of Liparis with three dorsal radials with notches and a ventral radial (Orr & Busby 
2006), as is true also for Prognatholiparis ptychomandibularis (Fig. 13H). Acantholiparis opercularis (Fig. 131) 
and all members of the clade Elassodisca share a simple, derived morphology of four large rounded radials, lack¬ 
ing notches and fenestrae, equally spaced in the basipterygium in all except C. notosaikaiensis; the morphology is 
unknown in C. cyclocephalus. The girdle is further derived in the clade Osteocareprocta with members having only 
two round dorsal radials (Kido 1988; Andriashev & Stein 1998; Duhamel et al. 2010). In the subgenus Temnocora, 
some differences are evident between the species: Careproctus comus, C. faunus, and C. candidus share the con¬ 
figuration of three dorsal and one ventral radials, one or two slightly notched with only two small fenestrae (Orr 
& Maslenikov 2007), while C. staufferi has the same number of radials that are strongly notched with three large 
fenestrae (Orr 2016). The derived Paraliparis penicillus has four small round radials in a basipterygium lacking 
fenestrae (Baldwin & Orr 2010), a common pattern found among members within the clade Paraliparia. 

In the clade Paracareprocta, pectoral girdle morphology is known for 14 of the 21 species. Lopholiparis flerxi 
is placed as the sister of other members of clade Paracareprocta in the COI tree and as the sister of other members 
of the subgenus Temnocora in the RADseq tree. Although its specific position is uncertain, its pectoral girdle is 
strikingly similar to the plesiomorphic condition of the girdles of Eumicrotremus and Liparis, having large radi¬ 
als with large fenestrae (Orr unpublished, based on a disarticulated specimen, UW 113885). Among species of the 
subgenus Careproctus, all species have three radials in the dorsal half of the basipterygium widely separated from 
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the ventral fourth, and at least two and usually three fenestrae (Andriashev 1986; Kido 1988; Chernova 2015; Orr 
et al. 2015). All species of the Melanura for which the pectoral girdle morphology is known— C. melanurus, C. col- 
letti, C. furcellus, C. cypselurus, C. micropus , and C. kidoi —have a similar pattern of three dorsal and one ventral 
radials (Chernova 2005b; Knudsen & Moller 2008; Maslenikov et al. 2015; Fig. 13J). Each radial is round and the 
girdle either lacks all fenestrae, as in C. micropus and C. kidoi (Knudsen & Moller 2008), or has a single fenestra 
between the scapula and first proximal radial. Paraliparis dactvlosus also has 3+1 round radials lacking notches and 
fenestrae except for the most dorsal (Fig. 13K), similar to species of Melanura. 

Tooth patterns. Stein (2012) discussed the utility of morphological characters among species of Paraliparis 
and identified tooth patterns, particularly uniserial simple teeth, as possibly having phylogenetic utility. The plesio- 
morphic condition among liparids is simple teeth (Knudsen et al. 2007), as found in cyclopterids and in Nectoliparis 
(Table 3). Among species of Liparis, nearly all species have trilobed teeth; the exceptions being the species formerly 
placed in the genus Polypera, P. greeni and P. shimushirae, which have mainly simple teeth, some with small lateral 
lobes (Burke 1930). Outside of Liparis, simple teeth are present in all members of the Pseudoliparia, some members 
of the Bathyphasmia, all members of the Paraliparia, and several members of the Aenigmoliparia (Table 3). Among 
clades within the Aenigmoliparia, trilobed teeth characterize all members of subgenus Temnocora and most mem¬ 
bers of clade Elassodisca, as well as most members of Allocareproctus and species placed incertae sedis (Table 3). 
Likely unique among Paraliparis, P. dactylosus has trilobed teeth (Table 3). 







FIGURE 13 A-G. Pectoral girdles of selected species of the Cyclopteridae and Liparidae: A) Eumicrotremus orbis, UW 1 11284; 
B) Nectoliparis pelagicus, UW 119455; C) Liparis gibbus, UW 119092; D) Crystallichthys cyclospilus, UW 47840; E) Care- 
proctus macrodiscus, FAKU 137835; F) Careproctus marginatus, FAKU 144616; G) Careproctus roseofuscus, FAKU 144615 
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In most liparids, teeth are aligned in narrow to broad oblique rows on both upper and lower jaws (Burke 1930; 
Table 3) and are uncommonly in a single series. Of the species Stein (2012, his table 4) identified as having uniserial 
teeth in the Ross Sea, all were placed within the subgenus Paraliparis and most within its derived clade as sister of 
P. bathybius, a North Atlantic species with uniserial teeth. Outside of this clade, species with uniserial teeth in our 
material are also found in the eastern North Pacific (P. rosaceus and P. paucidens: Stein 1978; Busby & Cartwright 
2009; Orr et al. 2005), Atlantic (P copei), and southeast Australia (Stein et al. 2001). Of these, P. rosaceus and P. 
copei are sister species within the P. copei species group, accompanied by three southern species, all presumed to 
have uniserial teeth and identified as P. neelovi (Duhamel et al. 2010), Paraliparis sp. from Australian Antarctica 
(FOAG559), and P. “bathybius ’’ (Rock et al. 2008). Paraliparis bathybius is known only from the Arctic and North 
Atlantic and has numerous small teeth (Andriashev 1954), and many specimens of R bathybius from the North At¬ 
lantic were recovered as the sister group of Paraliparis from the Ross Sea with uniserial teeth (Fig. 11, Clade M). 
In contrast, P. “bathybius” from the Antarctic was included in a clade of species with uniserial teeth from the North 
Atlantic, North Pacific, and Antarctic (P. copei, P. rosaceus, and P. neelovi, respectively; Fig. 11, Clade J). 



FIGURE 13 H-O. Pectoral girdles of selected species of the Cyclopteridae and Liparidae: H) Prognatholiparis ptychomandib- 
ularis, UW 156749; I) Acantholiparis opercularis, UW 118624; J) Careproctus sp. cf. melanurus, UW 119240; Yd) Paraliparis 
dactylosus, UW 116232; L) Rhinoliparis attenuatus, UW 113736; M) P cephalus, UW 117527; N) Paraliparis pectoralis, UW 
117515; O) P. ulochir, UW 150802. 

Additionally, a species of the eastern North Pacific identified as Paraliparis sp. (Steinke et al. 2009) with simi¬ 
lar morphology to P. rosaceus is clearly distinct in COI data (Steinke et al. 2009). It was placed within the subgenus 
Paraliparis, as sister of the southern species P voroninorum and P. neelovi and near P. paucidens, all with uniserial 
teeth. Exceptionally within this clade is P. mawsoni, a species with the plesiomorphic condition of tooth bands (An¬ 
driashev 1986), placed uncertainly in the polytomy near the base of the clade (Fig. 11, clade M). 
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Postorbital cephalic pore. Chernova (2005b) noted for the subgenus Careproctus the absence of the postor¬ 
bital pore in North Atlantic and Arctic species and its presence in Southern Hemisphere species. Among members of 
the subgenus Careproctus in our analysis, Orr et al. (2015) suggested that two clades would form the Careproctus 
rastrinus complex of North Pacific and Arctic species based on the presence or absence of the postorbital pore and 
suggested that C. acanthodes may be more distantly related. Our results support the presence of two clades, one with 
C. phasma and C. lerikimae closely related to C. reinhardti and all lacking the postorbital pore; the other with the 
five remaining species, including C. acanthodes, with postorbital pore present. Members of the third major clade of 
the subgenus, C. gilberti and C. ostentum, have a postorbital pore. 


Zoogeography and ecology 

Distributions of species in trees. The sister groups of all other liparids are species from the North Pacific, including 
the highly apomorphic Nectoliparis and most species of the plesiomorphic Liparis, both genera that are limited to 
the Northern Hemisphere. Among species of Liparis, North Pacific species form the sister groups of North Atlantic 
species. The North Atlantic species L. atlanticus, L. inquilinus, L. liparis, and L. montagui were recovered in a de¬ 
rived clade among the Arctic species L. fabricii, L. tunicatus, and L. bristolensis. Among western Pacific species, 
L. curilensis was placed within a clade of the eastern North Pacific species L. callyodon, L. florae, and L. greeni', 
L. punctulatus with L. rutterv, and L. tessellatus with L. pulchellus, while a clade of L. agassizii, L. tanakae, and L. 
chefuensis was the sister of a poorly supported clade of the eastern North Pacific and Arctic species L. dennyi, L. 
gibbus, and/,. bathyarcticus. 

Although the relationships among more derived clades are uncertain, the mix of zoogeographic regions repre¬ 
sented within most clades indicates multiple invasions of the Southern Ocean from the South Pacific, and the North 
Atlantic from the Southern Ocean and the North Pacific. Southern species of Careproctus were found within four 
clades placed in several positions among North Pacific or deep-water species. 

Among species of the deep-water clade Pseudoliparia, the southern abyssal species C. crozetensis is the sister of 
the deep abyssal and hadal genera Notoliparis and Pseudoliparis from the Mariana, Kermadec, and Japan trenches. 
Species of Notoliparis are known from the southern Kermadec, Macquarie-Hjort, South Orkney, and Peru-Chile 
trenches (Andriashev 1975, 1978; Stein, 2005; Stein 2016; Linley et al. 2016); Pseudoliparis had been known only 
from the northern Japan and Kuril-Kamchatka trenches (Andriashev 1955, Andriashev & Pitruk 1993) until the 
recent description of P. swirei from the southern Mariana Trench (Gerringer et al. 2017a). 

fn the Bathyphasmia, the southern species C. georgianus was sister of the eastern North Pacific C. ovigerus. 
Within the large polytomy of the clade Aenigmoliparia, three southern species—C. continentalis, C. discoveryae, 
and C. longipectoralis —formed a clade (subgenus Careproctula ) with the northern Osteodiscus cascadiae. A clade 
of one or two southern species (both identified as Careproctus sp.) was sister of an undescribed northern species 
(Orr, unpublished) within the subgenus Allochir. 

The Crystallichthyia, a clade of northern species of Careproctus and Crystallichthys cyclospilus that are all 
exclusively from the western North Pacific except for Crystallichthys, was the sister group of the clade Paraliparia, 
the large clade of species and genera that lack a pelvic disc. Although its position was poorly supported in the COl 
tree, in the RADseq tree it is highly supported. Thus, the derivation of the disc-less clade from the North Pacific is 
likely. 

In contrast with the scattered southern species of Careproctus, all southern species of Paraliparis were placed 
within the Paraliparia. Within this clade, southern species were spread over at least four internal clades, plus the 
single P. thalassobathyalis. Three of these clades also included both North Atlantic and North Pacific species; one 
included a single North Atlantic species. Three other northern species of Paraliparis — P. dactylosus, P. grandis, and 
P. penicillus —were placed among northern species of Careproctus within the Aenigmoliparia. 

Zoogeography. The Liparidae likely arose in the North Pacific in the late Eocene-early Miocene about 45 to 
25 mya (Stein et al. 2001; Near et al. 2012; results of Shen et al. 2017 who estimated 28-10 mya were inappro¬ 
priately calibrated by a fossil of the derived recent species Gasterosteus aculeatus from 13.3 to 13.0 mya), with a 
burst of rapid speciation that began at about 30 mya (Rabosky et al. 2018). Other than a single questionable fossil 
from the Oligo-Miocene of Belgium (Chernova et al. 2004), fossils of liparids are unknown (Andriashev & Stein 
1998). Our results indicate that liparids were derived from a shallow-water ancestor in the North Pacific, likely the 
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eastern North Pacific. The highly apomorphic pelagic species Nectoliparis pelagicus of the North Pacific diverged 
early. One, possibly two, lineages of Liparis diversified into shallow waters. Following this split, a lineage diverged 
into deep waters and from this deep-water group all other liparids likely evolved, based on the deep-water group’s 
sister group position to most other liparids and morphological evidence from the plesiomorphic pectoral girdle, 
among other characters (Andriashev 1978, 1998; Stein 2016). Deep-water groups further diversified throughout the 
northern and southern Pacific into the relatively cooler shallow waters at the poles, persisting and further diversify¬ 
ing in colder waters as the Miocene climate generally warmed 17-15 mya. Beginning no earlier than about 41 mya, 
via the opening of the Drake Passage in the Miocene about 41-20 mya (Scher & Martin 2006), liparids dispersed 
into the South Atlantic in multiple invasions through the Southern Ocean and around Antarctica (Andriashev 1986; 
Andriashev & Stein 1998; Stein et al. 2001). Eastern South Pacific species may also have dispersed into Australian 
waters (Stein et al. 2001). 

Clades from both the North Pacific and Southern Hemisphere likely invaded the North Atlantic more recently. 
Paraliparis and perhaps Careproctus dispersed north out of the South Atlantic and Southern Ocean (Andriashev 
1991, 1997; Andriashev & Stein 1998; Chernova 1999), similar to “ambiguous” patterns evident within the mac- 
rourid genus Corvphaenoides (Gaither et al. 2016). The Antarctic Intermediate Water Mass may have provided a 
mechanism for this dispersal north (Koslow et al. 1994; Stein et al. 2001; Probert 2017). Eastern North Pacific 
species are closely related to Antarctic species in some clades (e.g., Bathyphasmia, Allochir, subgenus Paraliparis) 
and members of these clades of bathyal species may have also dispersed across the Central American Seaway over 
the open Panamanian isthmus (Balushkin 2012), where at its southern extent a complete interchange of Atlantic and 
Pacific waters at bathyal depths was occurring before 20 mya until 12 mya at the early collision of the American 
plates (Coates et al. 2004). This region shallowed to 200 m or less at 11-10 mya (Sepulchre et al. 2014; Montes et 
al. 2015) before its complete closure in the Pliocene 2.76-2.54 mya (Sepulchre et al. 2014). Other interpretations 
visualized a dispersal of North Pacific species along the coast of South America and into the Antarctic and South 
Atlantic (Andriashev 1987, 1993; Balushkin & Voskoboinikova 2008), as well as to Australia (Stein et al. 2001). 
Examination of molecular and other data of uncommon Pacific South American species will be important to better 
understand these distributional patterns, especially their relationship to Australian species (Stein et al. 2001). 

To the north, shallow-water species of Liparis and Careproctus were also able to invade the North Atlantic since 
about 7.4-4.8 mya in the upper Miocene, and perhaps earlier in the Middle Miocene during two periods 17-16 mya 
and 12—11 mya (Sher 1999), across the open but shallow Bering Strait, which has likely never been deeper than 
50 m (Marincovich and Gladenkov 2001). This shallow water excluded deep-water liparid genera such as Para¬ 
liparis from dispersing to the Arctic and North Atlantic via this route (Andriashev 1986, 1991). Movements from 
the North Pacific into the Arctic and North Atlantic are well-established in the mollusk fossil record (Marincovich 
& Gladenkov 1999) and through molecular and phylogenetic analyses of many other fish families (e.g., Cottidae: 
Knope 2013; Scorpaenidae, Sebastinae: Hyde & Vetter 2007; Zoarcidae: Anderson 1994; Moller & Gravlund 2003), 
although some families, including the cyclopterids (Ueno 1970) and gadids (Coulson et al. 2006; Carr & Marshall 
2008; Roa-Varon & Orti 2009), may have originated in the Atlantic and diverged into the Pacific. 

Reproductive parasitism. All liparid eggs identified to species taken from inside crabs in the Northern Hemi¬ 
sphere have been of species in the Clade Melanura. Gardner et al. (2016) previously identified eggs of Careproctus 
melanurus, C. colletti, C. furcellus, and C. simus taken from Golden and Scarlet King Crabs (Lithodidae: Lithodes 
aequispinus and L. couesi). Two sequences from GenBank were from eggs deposited in Blue King Crabs ( Paral- 
ithodes platypus ; Yanagimoto 2015) and based on our data could be identified as eggs of species of the subgenus 
Allochir. one sequence (LC002654.1) was identical to sequence data of C. furcellus-, the other (LC002656.1) was 
identical to sequence data of C. colletti. Although Yanagimoto (2015) conducted a GenBank BLAST analysis to 
identify the eggs, without success, and discounted C. furcellus as an identification, sequences labeled as C. furcellus 
in GenBank at the time were based on misidentified specimens of C. cypselurus (Appendix Tables 1 and 2), and 
data for C. colletti were not yet publicly available. Yanagimoto’s (2015) report is the only published documentation 
of reproductive parasitism in Blue King Crabs. Other instances of parasitism have been reported from the Southern 
Hemisphere (citations in Gardner et al. 2016) but the only southern species known to deposit eggs in crabs are C. 
albescens (= C. griseldea ) and C.falklandicus (Balbontin et al. 1979; Melville-Smith & Louw 1987; Andriashev & 
Prirodina 1990; Chernova 2014) and neither were available for our analysis. 

In contrast to the eggs of Liparis, which are commonly found in shallow waters, identified egg masses of 
deep-water snailfishes are rare and few have been reported outside of crabs (e.g., C. rastrinus: Chernova 1992; Al- 
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locareproctus unangas : Busby et al. 2006; Squaloliparis dentatus : Poltev & Stelcsova 2010; C. fulvus: Chernova 
2014; Elassodiscus tremebundus, Paraliparis rosaceus: Overdick et al. 2014; Gardner et al 2016). Except for eggs 
identified as Careproctus sp. (as either C. colletti or C. melanurus) that were found in trawls with lithodid crabs 
(Overdick et al. 2014), all species depositing eggs outside crabs, and for which we had molecular data, were out¬ 
side the Melanura clade. Careproctus fulvus, for which we did not have molecular data and which deposits eggs in 
sponges, is a member of the subgenus Careproctus (Chernova 2014). 


Conclusions 

Our phylogenetic analysis of nearly a third of the species and over half the genera of the Liparidae from both South¬ 
ern and Northern hemispheres has provided useful insights to continue to build on our understanding of the evolu¬ 
tion of the family. However, the relationships of this diverse family of more than 400 species and 32 genera remain 
in large part unclear. Our results from a single-gene analysis of a wide range of taxa, bolstered by the data-rich and 
largely congruent genome-wide analysis of North Pacific species, have provided confidence in many nodes of the 
phylogeny, but several important transitions in liparid evolution remain unresolved. 


TABLE 4. Species of Liparidae and outgroups used to construct a RADseq phylogenetic tree, listed as clades and 
species appear in the tree. Letter designations of clades and taxonomic names are referenced in the text and in Fig¬ 
ures 8-12. Species in bold were placed in an alternative position in the COI tree. 


Lettered clade 

Rank 

Species 

Position in COI tree 


Cyclopteridae 

Eumicrotremus orbis 



Liparidae 
genus Liparis 

Nectoliparis pelagicus 


Clade Q 

subgenus Polypera 

Liparis florae 



subgenus Careliparis 

L. gibbus 



Liparis rutteri species group 

L. rutteri 


Clade D 

Clade Pazlliparia 



Clade E 

Clade Aenigmoliparia 

Prognatholiparis ptychomandibularis 

Aenigmoliparia 
incertae sedis 


Clade Elassodisca 

Paraliparis penicillus 

Aenigmoliparia 

incertae sedis 



Elassodiscus caudatus 



unnamed clade 

Aliocareproctus spp. (A. jordani, A. ungak) 

Aenigmoliparia 

incertae sedis 


subgenus Temnocora 

Lopltoiiparis fierxi 

Careproctus comus 

Careproctus faunus 

Careproctus staufferi 

Careproctus candidus 

Paracareprocta 

. continued on the next page 
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TABLE 4. (Continued) 


Lettered clade 

Rank 

Species 

Position in COI tree 

Clade G 

Clade Paracareprocta 



Clade P 

P. dactylosus species group 

Paraliparis dactylosus 


Clade 0 

subgenus Careproctus 

C. gilberti species group 

Careproctus gilberti 



C. rastrinus species group 

Careproctus scottae 



C. reinhardti species group 

Careproctus phasma 


Clade N 

Clade Melanura 




subgenus Allochir 

Careproctus cvpselurus 

Careproctus melanurus 

Careproctus sp. cf. melanurus 
Careproctus furcellus 


Clade F 

Clade Ambidisca 



Clade H 

Clade Crystallichthyia 

Crystallichthys cyclospilus 


Clade Q 

Clade Paraliparia 

Lipariscus nanus 

Rhinoliparis attenuatus 


Clade 1 

P. ulochir species group 

Paraliparis ulochir 



Although COI has become a widely accepted means for identifying species and describing differences among 
them, RADseq data contains an order of magnitude more sequence data. It is not unprecedented for RADseq trees 
and COI trees to differ. In some cases, RADseq trees have been considered an improvement (Razkin et al. 2016), the 
RADseq data providing resolved trees for cases that mtDNA (and nuclear DNA) could not (Diaz-Arce et al. 2016), 
as in the case with liparids. While other work has concluded that reconstruction of recently diverged lineages may 
still be difficult with RADseq data in some cases (Suchan et al. 2017), this is unlikely to be the case for the more 
distantly related species of our liparid dataset. 

While our analysis of COI sequence data successfully resolved many shallow and deep nodes of the global li¬ 
parid tree, it may be an inadequate source of data for this group because of multiple superimposed nucleotide substi¬ 
tutions in deep phytogenies. Our ability to select loci containing an appropriate level of variation with the RADseq 
analysis provided a completely resolved tree with high support for a much smaller group of species from the eastern 
North Pacific—expanding taxon sampling for additional RADseq analyses will be an important avenue for further 
research. These molecular phylogenies provide a skeleton of proposed relationships that will allow detailed phylo- 
genetically based morphological assessments, leading to taxonomic and nomenclatural clarity. 
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Careproctus colletti UW 118684* MH630284 SGJWO049-18 

Careproctus colletti UW 118687* MH630281 SGJWO053-18 

Careproctus colletti UW 119297* MH630280 SGJW0050-18 

Careproctus colletti UW 150595 JQ354031 FMV577-11 Careproctus sp. 

Careproctus comus UW 119197* MH630285 SGJW0060-18 
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Careproctus marginatus FAKU 144616 LC337285 N/A 

Careproctus melanurus UW 115145 JQ354030 FMV419-09 

Careproctus melanurus UW 150589* MFI630295 SGJWO048-18 

Careproctus sp. cf. melanurus UW 118920* MFI630296 SGJWO047-18 

Careproctus micropus TMUAL-9299 N/A DSFNG033-11 
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Liparis dennyi UW 48801/SIO 05-195 JQ354178 FMV134-08/MFC353-08 

Liparis dennyi UW 155675* MH630319 SGJWO014-18 

Liparis dennyi UW 119994 JQ354188 FMV529-11 Liparis sp. 

Liparis fabricii ZMUB Vp2013-66 N/A DSFIB528-14 

Liparis fabricii ZMUB 10313 N/A DSFIB751-16 
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Liparis pulchellus RBCM NEOCAL07-0036 FJ164724 TZFPA183-07 

Liparis pulchellus UW 112760 JQ354184 FMV249-08 

Liparis pulchellus RBCM NEOCAL07-0039 FJ164721 TZFPA186-07 

Liparis pulchellus UW 110220/SIO 05-200 GU440378 FMV218-08/MFC354-08 

Liparis pulchellus RBCM NEOCAL07-0032 FJ164725 TZFPA179-07 
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Nectoliparis pelagicus RBCM NEOCAL07-0046 FJ164907 TZFPA193-07 

Nectoliparis pelagicus RBCM NEOCAL07-0047 FJ 164906 TZFPA 194-07 

Nectoliparis pelagicus UW 117244-03* MH630323 SGJWO018-18 

Notoliparis kermadecensis HADES 100340 KY659176 N/A 

Notoliparis steward HADES 100344 KY659177 N/A 
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Careproctus furcellus UW 117984* MH715524 N/A 

Careproctus furcellus UW 118636* MH715525 N/A 

Careproctus furcellus UW 118724* MH715526 N/A 

Careproctus furcellus UW 119295* MH715527 N/A 

Species Catalog number Genbank accession BOLD sequence Original identification 
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Careproctus gilberti UW 152425-01* MH715552 N/A 

Careproctus gilberti UW 153160-01* MH715553 N/A 

Careproctus gilberti UW 153160-02* MH715554 N/A 

Species Catalog number Genbank accession BOLD sequence Original identification 
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Careproctus reinhardti TMU AL-9252 N/A DSFNG022-11 

Careproctus reinhardti TMU AL-9254 N/A DSFNG023-11 
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Species Catalog number Genbank accession BOLD sequence Original identification 








3 

bo 

a 

ex 

8 

3 

-Si 


3 

5 

S' 

<3 


- 

ON 

<n 




'sj- 

■rj- 

'3- 

NO 

cn 

1 

o 








f-* 

1—1 

o 

CN 

'sf 

'sf 

'sf 

1 

i 

c-1 

ON 

<N 





in 

cn 


o 

nO 

ON 




<N 

cn 

cn 


o 

z 

Uh 

o 

r- 

o 

(N 

<n 

in 

<n 

7 

— ] 

ss 

R 

2 

r- 

o 

R 

NO 

R 

NO 

R 

2 

E 

s 

E 

s 

E 

§ 

(/) 

C/} 

CQ 

R 

R 

R 

c/3 

GO 

GO 

CQ 

Q 

Q 

2: 

o 

o 

o 

Q 

Q 

Q 

O 


Z 


z 


z 


r- 

ON 


< 

Z 


z 


z 


oo 

oo 

ON 

CN 


a 


< 

z 


z 


oo 

© 

r- 

NO 

00 

QQ 

CLh 

R 

N 

H 


z 


OO 

NO 

00 

f-* 

ON 

o 

—t 

CN 

cn 

Tf 

N- 

-'t 

"sf 

<n 

in 

m 

cn 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

CN 

cn 

cn 

cn 

cn 

cn 

cn 

cn 

i —i 

<n 

in 

m 

in 

in 

in 

in 

o ^ 

i? 5: 5: 5: i? p 

o 

X 

o 

E 

o 

D 

o 

X 

o 

P 

o 

X 




'si- 

r-* 

m 

in 

t> 

X 


S 


«n no 
r- r- 


r- r- 

X ffi 


r- oo 

r- r- 


r- r- 

X ffi 


ON 

r- 


r- 

X 


© 


00 00 00 
m ^ «n »n m 
in m »n 


cn 

in 

o 

X 

*4 


t— t— 

X ffi X 


ON 

n 

"sf 

■'T 

NO 


CN in CN 
Tj- h °o 





C"- 

f~-. 

o 

NO 

CN 

NO 

cn 

cn 

cn 


CN 

cn 

NO 

00 


CN 

* 







* 

* 

CN 

* 

cn 


00 

CN 


«n 

ON 

o 

R 

Oh 

Oh 

o 

CN 

o 

CN 

o 

CN 

NO 

© 

NO 

© 

NO 

O 

NO 

o 

NO 

© 

CN 

© 

CN 

o 

in 

* 

* 

r- 

* 

00 

* 

OO 

* 

o 

o 

CN 

O 

CN 

o 

CN 

00 

o 

ON 

o 

L 

R 

J 

R 

CX 

Q. 

CX 

CN 

O 

o 

O 

o 



cn 

cn 

cn 

cn 

r- 

ON 

N" 

00 

00 

00 


J 

o 

Cn 

CN 

O 

O 

O 

> 

> 

> 

in 

ON 

ON 

ON 

ON 

On 

ON 

ON 

o 

'sf 

N - 

r- 

r- 

’-H 

N - 

's|- 

's|- 

_' 

AJ 

00 


cn 

cn 

cn 

cn 

cn 

r-- 

r- 

l> 

NO 

ON 

ON 

On 

ON 

cn 

N" 

N" 

N" 


< 

's|- 

03 

00 

03 

03 

03 

03 

03 

*—i 

l-H 

7 —1 

l-H 

7—* 

7 —1 

7 —1 

7 —1 

—H 


7 —1 



in 

i n 

«n 

IT) 


1 —) 

m 

R 



R 

R 

R 


i-H 

’“H 

i-H 

’“H 

l-H 

’“H 

i-H 

’“H 

i-H 

’-H 

’-H 

’-H 

i-H 

’-H 

’-H 

’-H 

H 


§ 

z 

i— 1 

s 

R 

s 

§ 

i—' 

§ 

R 

s 

R 

§ 

r 

§ 

£ 

£ 

£ 

£ 

£ 

£ 








£ 




U 

CQ 

H 

N 

N 

N 

N 

N 

N 

N 

N 

P 

P 

P 

P 

P 

p 

P 

p 

P 

P 

P 

P 

P 

P 

P 

P 

P 



-d 

a i 


o 


U 


w 


is 

is 

is 

is 

is 

is 

is 

is 



















p 

ea 

Q 

Q 

Q 

Q 

Q 

Q 

s. 

Q 

s. 

Q 

s. 

Q 

§ 

§ 

§ 

§ 

§ 

§ 

§ 

SO 

§ 

03 

3 

03 

3 

03 

3 

03 

3 

03 

3 

03 

3 

03 

3 

03 

3 

03 

3 

d. 


'S 

'g 

s; 

'g 

'g 

’p 

Si 

'g 

.s; 

’p 

Si 

'g 

Si 

'g 

0 

O 

O 

0 

O 

O 

O 

O 

•1 

•§ 

•1 

•1 

•1 

•1 

•1 

•1 

g 

C /3 

H 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 


03 

03 

03 

03 

03 

03 

d 

X 

03 

03 

to 

to 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 

03 


S 

S 

= 

a 

s 


s 


S 

S 

S 

3 


3 

.3 

3 

,3 

3 

3 

3 


3 


3 

3 

-g 

a 

Q 


CJ 

U 

o 

NO 

NJ 

1 ) 

NJ 

NJ 


NJ 

1 ) 

NJ 

13 

NJ 

13 

SO 

13 

3 

13 

3 

3 

3 

3 

3 

3 

3 

Z 

g 

P 

g 

g 

g 

g 

g 

g 

g 

g 

P 

g 

g 

g 

P 

g 

g 

g 

P 

g 

g 

g 

g 

g 

g 

g 

g 

W 


5 , 


a. 






sx, 

5 , 




5 , 



EX, 

5 , 



EX, 

EX, 

EX, 

EX, 

EX, 

EX, 

CLh 

g 

g 

g 

g 

g 

g 

g 


g 

g 




g 

g 

g 

g 

g 

g 

g 

g 

g 

g 

g 

g 

g 

g 

Pm 

Q 


■3 

a 



Q 



Q 

Q 

Q 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


U 

o 

o 

O 

o 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 


MOLECULAR PHYLOGENETICS OF THE LIPARIDAE 


Zootaxa 4642 (1) © 2019 Magnolia Press • 67 


Careproctus n. sp. RBCMINV792 FJ164430 TZFPA034-06 Careproctus canus 

Careproctus n. sp. RBCM 15946 FJ164431 TZFPA070-06 Careproctus canus 

Crystallichthys cyclospilus UW 48076* MF1715583 N/A 

Species Catalog number Genbank accession BOLD sequence Original identification 
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Liparis dennyi UW 119996 JQ354190 FMV531-11 Liparis sp. 

Liparis dennyi UW 155676* MH715616 N/A 

Liparisfabricii ARC 26710 KCO15559 SCAFB242-07 

Species Catalog number Genbank accession BOLD sequence Original identification 
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Liparis gibbus UW119134 N/A DSFIB068-11 

Liparis gibbus ZIN 54878 HM421748 DSFAL608-09 

Liparis liparis MT 1807 KJ204968 BNSF316-11 

Species Catalog number Genbank accession BOLD sequence Original identification 
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0RR ET AL. 


Nectoliparis pelagicus RBCM NEOCAL07-0045 FJ164908 TZFPA192-07 

Nectoliparis pelagicus UAM PSR 2143-01 N/A DSFIB188-12 

Nectoliparis pelagicus UAM PSR 2143-02 N/A DSFIB189-12 

Species Catalog number Genbank accession BOLD sequence Original identification 
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APPENDIX TABLE 3. Material examined in 31 groups of species designated for RADseq analysis. SRA = National 


Center for Biotechnology Information sequence read archive. 

Group 

Catalog number 

Species 

SRA sample number 

1 

UW 119822 

Eumicrotremus orbis 

Eorbis822 

2 

UW 119829 

Lopholiparis jlerxi 

Lflerx9829 

3 

UW 119776 

Paraliparis ulochir 

Pulochir76 

4 

UW 156086-1 

Allocareproctus jordani 

Ajordani86 

4 

UW 156089 

Allocareproctus jordani 

Ajordani89 

4 

UW 155938 

Allocareproctus ungak 

Ajordan938 

5 

UW 112756 

Liparis rutteri 

Lrutteri56 

6 

UW 49350-2 

Elassodiscus caudatus 

Ecaudatus2 

6 

UW 119593-1 

Elassodiscus sp. cf. caudatus 

Ecaudatus3 

6 

UW 155734 

Elassodiscus sp. cf. caudatus 

Ecaudatus4 

6 

UW 49350-1 

Elassodiscus caudatus 

Ecaudatusl 

7 

UW 152008 

Paraliparis dactvlosus 

Pdactyl08 

7 

UW 154451 

Paraliparis dactvlosus 

Pdactyl51 

8 

UW156091-1 

Careproctus candidus 

Ccandidus 

8 

UW 156091 

Careproctus candidus 

Ccandidu2 

9 

UW 116030 

Careproctus cypselurus 

Ccypselur 

10 

UW 156087-2 

Careproctus sp. cf. melanurus 

Ccfmelanu2 

11 

UW 152476 

Lipariscus nanus 

Lnanus2476 

12 

UW 119187 

Crystallichthys cyclospilus 

Crystall87 

12 

UW 48076-2 

Crystallichthys cyclospilus 

Crystalli2 

12 

UW 48076-1 

Crystallichthys cyclospilus 

Crystallil 

13 

UW 119291 

Careproctus bowersianus 

Cbowersian 

14 

UW 151304 

Careproctus scottae 

Cscottae04 

14 

UW 151302 

Careproctus scottae 

Cscottae02 

15 

UW 155711 

Careproctus staufferi 

Cstauffel 1 

15 

UW 119196 

Careproctus staufferi 

Cstauffe96 

16 

UW 119197 

Careproctus comus 

Ccomus970 

16 

UW 119197-1 

Careproctus comus 

Ccomus971 

17 

UW 119821 

Eumicrotremus orbis 

Eorbis821 

17 

UW 151223 

Eumicrotremus orbis 

Eorbis223 

17 

UW 152461 

Eumicrotremus orbis 

Eorbis461 

18 

UW 151666 

Liparis florae 

Lflorae666 

19 

UW 116036 

Prognatholiparis ptychomandibularis 

Pptychoman 

20 

UW 154481-1 

Careproctus gilberti 

Cgilberti81 

21 

UW 117244 

Nectoliparis pelagicus 

Npelagicus 

22 

UW 119176 

Liparis gibbus 

Lgibbusl76 

22 

UW 119043 

Liparis gibbus 

Lgibbus043 

23 

UW 156084 

Careproctus faunus 

Cfaunus841 

24 

UW 113696 

Allocareproctus ungak 

Aungak696 

25 

UW 119298 

Careproctus furcellus 

Cfurcellu98 

26 

UW 151307 

Careproctus phasma 

Cphasma07 

26 

UW 151261 

Careproctus phasma 

Cphasma61 

27 

UW 119192 

Paraliparis penicillus 

Ppenicillu 

28 

UW 115871 

Rhinoliparis attenuatus 

Rattenuatu 

29 

UW 150589 

Careproctus melanurus 

Cmelanu89 

30 

UW 150588 

Careproctus melanurus 

Cmelanu88 

31 

UW 118920 

Careproctus sp. cf. melanurus 

Ccfmelanul 
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